U.S. DEPARTMENT OF COMMERCE 
National Technical Information Service 


I 

i 

'6 

j 

i 

k ! 

I 

J 

rl 

>1 


N77-23211 i 

I 

ANALYSIS OF SHEAR TEST METHOD FOR COMPOSITE 
LAMINATES 

HENRY W, BERGNER, ET AL e 

VIRGINIA POLYTECHNIC INSTITUTE AND STATE ? 

UNIVERSITY 

BLACKSBURG> VIRGINIA r 



APRIL 1977 


2 


PlSt'mgT^r no^[ 

Appioved to; puoiic tsi«<u«( 1 
Piimbunon Didaiutad I 



DTIC QUALITY INSPECffiDS 



^H17- dd£// 

College of Engineering 

Virginia Polytechnic Institute and State University 
Blacksburg, VA. 24061 


VPI-E-77-14 


April, 1977 


Analysis of Shear Test Method 
for Composite Laminates 

Henry W. Bergner, JrJ 

2 

John G. Davis, Or. 

3 

Carl T. Herakovich 


Department of Engineering Science and Mechanics 
Interim Report Number 6 

NASA-VPI&SU Cooperative Program in COTiposite 
Materials Research and Education 


Accesion For . 

NTiS CRA&I B 

DTiC TAB □ 

U nan r>ou need □ 

Justification 


By 

Distribution / 


Avaiiabiiiiy Codes 


Dist 

Avaii and/or 

Special 

A'l 



NASA Grant NGR 47-004-129 

Prepared for: Materials Application Branch 

National Aeronautics & Sapee Administration 
Langley Research Center 
Hampton, VA. 23665 

"* Graduate Student 
2 

Materials Engineer - NASA-Langley 

3 

Professor 


Approved for public release, distribution unlimited. 


REPRODUCED BY 

NATIONAL TECHNICAL 
INFORAAATION SERVICE 

U, S. DEPARTMENT OF COMMERCE 



i 


BIBLIOGRAPHIC DATA 1. Report No. 1. 

sheet VPI-E-77-14 

3, Recipient's Accession No. 


4. Title and Subtitle 

1 ANALYSIS OF SHEAR TEST METHOD FOR COMPOSITE LAMINATES 

! 


6. 

7. Author(s) 

Henry W. Beroner. Jr., John G. Davis, Jr. Carl T. Herakovic 

8. Performing Organization Rept. 

1 VPI-E-77-14 

!?. Performing Organization Name and Address 

1 Virginia Polytechnic Institute and State University 

Engineering Science and Mechanics 
1 Blacksburg, Virginia 24061 

10. Ptoject/Task/Work Unit No. 

11. Contract/Grant No. 

NASA NGR-47-004-129 

iiJ, Sponsoring Organization Name and Address 

i National Aeronautics & Space Administration 

1 Langley Research Center 

i Hampton, Virginia 23665 

1 

13. Type of Report & Period 
Covered 

14. 


Supplementary Notes 


f 16. Abstracts 

I 

see page n 


17. .Key Words and Document Analysis. 17a. Descriptors 

Composites, In-plane Shear, Finite Elements, .Rail Shear, losipescu. Cross 
Beam, Slotted Coupon, Elastic Fixtures, Thermal Stresses, Graphite-epoxy, 
Graphite-polyimide 


?7b. Wentifiers/Open-Ended Terms 


! 17c. tOSATI Field/Group 


18. Availability Statement 

Distribution unlimited 

19. Security Class (This 
Report) 

UNCLASSIFIED 

21. No. of Pages 


20. Security Class (This 

“Unclassified 

22. Price 


FORM NTis-sa (REV. 10-731 ENDORSED BY ANSI AND UNESCO. THIS FORM MAY BE REPRODUCED uscomm.dc »2«5-P74 



















FOREWORD 


This report documents a portion of the work accomplished under NASA 
Grant NGR-47-004-129 during the period September, 1975 through April, 1977. 
Mr. Bergner was in residence at NASA-Langley during the period September 
1976 - February, 1977. 

The authors are grateful to Mr. Donald J. Baker of the U. S. Army Air 
Mobility Research and Development Laboratory at NASA's Langley Research 
Center, and Mr. Edward A. Humphreys and Ms. Frances Carter of Virginia 
Polytechnic Institute and State University for their contributions to this 


ANALYSIS OF SHEAR TEST METHODS FOR COMPOSITE LAMINATES 


ABSTRACT 


An elastic plane stress finite element analysis of the stress 
distributions in four flat test specimens for in-plane shear response 
of composite materials subjected to mechanical or thermal loads is 
presented. The shear test specimens investigated include: slotted 

coupon, cross beam, losipescu, and rail shear. Results are presented 
in the form of normalized stress contour plots for all three in-plane 
stress components. 

The slotted specimen is studied for three graphite-epoxy lami- 
nates ([0], [90], [±45]s); the cross beam is studied for two graphite 
epoxy laminates ([0/90]g, [±45]g) ; the losipescu and rail shear speci 
mens are studied with several materials (steel, graphite-epoxy, and 
graphite-polyimide) and several laminates ([0], [90], [0/90]5, [±45]j 
[0/90/±45]g) with rigid and elastic fixtures loaded mechanically or 
thermally. Geometric alterations are also investigated. 

The study shows that the cross beam, losipescu, and rail shear 
specimens have stress distributions which are more than adequate for 
determining linear shear behavior of composite materials. Laminate 
properties, core effects, and fixture configurations are among the 
factors which were found to influence the stress distributions. 
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2. LITERATURE REVIEW 


2.1 The Slotted Shear Specimen 

The slotted shear specimen has been used to measure shear proper- 
ties of both metals and composites. Fenn and Clapper [1] used the 
geometry of Fig. la for their studies on metals. Breindel et al [2] 
suggested a standardized metal specimen thickness, width, and length of 
shear path (Fig. la). losipescu's photoelastic study [3] of this 
specimen indicated that the state of stress for a slotted specimen is 
generally unacceptable for quantitative shear behavior due to the lack 
of pure uniform shear stress in the test area. 

Elkin et al [4] used the slotted specimen in Fig. lb for. shear 
tests on three unidirectional graphite-epoxies. They reported low 
strength values which were attributed to the stress concentrations at 
the edges of the holes. Their finite element analysis indicated a max1 
mum stress concentration of 1.57 at the edge of the holes. Other 
variations of the slotted shear specimen shown in Figs. Ic and Id have 
been used for composite materials [5]. 

2.2 The Cross Beam Shear Specimen 

The cross beam shear specimen has been used for determining linear 
and nonlinear shear stress-strain response as well as ultimate shear 
strength values for graphite-epoxy and boron-epoxy laminates. As shown 
in Fig. 2a, the cross beam shear specimen is loaded in positive and 
negative bending to produce a biaxial state of tension and compression 
over the test area at the center of the specimen. Stresses on 45° 
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FIGURE I. SLOTTED SHEAR SPECIMENS 
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planes are presumed to be pure, uniform shear over the test section. 
Laminates comnonly used with the cross beam are [0/90]^ and [±45]^. 

Petit [6] compared experimental results of this test method to 
those of the [±45]^ tension coupon and found good agreement for both 
linear and nonlinear shear stress-strain response of unidirectional 
laminates. Waddoups [7] compared experimental ultimate shear strength 
values obtained by this method to those of the short beam shear test and 
reported that the cross beam specimen produced considerably higher 
ultimate shear strength values than those of the short beam shear test. 

2.3 The losipescu Shear Specimen 

As mentioned previously, an extensive photoelastic study was con- 
ducted by losipescu [3] on some of the slotted shear specimens that had 
been used for shear testing. The results of this study generally 
indicated a need for a more precise shear test method. A new specimen 
(the losipescu specimen) was consequently developed as a test method for 
shear strength of metals. 

Although the losipescu shear specimen has not been used extensively 
for shear testing of composite materials, the encouraging results from 
metals indicate that this specimen should be analyzed for composite 
materials to evaluate any potential for use as a test method. The ef- 
fects of several modifications to the original losipescu specimen are 
studied for use with composite materials. While metal losipescu speci- 
mens are normally notched in two directions, the specimen in the present 
study is notched in only one direction (Fig. 2b) so that the laminate 
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FIGURE 2. VARIOUS SHEAR SPECIMENS 
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stiffness in the test section is uniform and strain measurements may be 
taken at the center. The specimen is loaded to produce a shear stress 
acting between the notch ends with zero bending moment and zero normal 
stress. The reduced specimen area between the notch ends ensures 
failure in this region. 

2.4 The Rail Shear Specimen 

The rail shear specimen has been used quite often for composites 
with various laminate configurations and at elevated temperatures as 
well as room temperature. As shown in (Fig. 3) and in studies by Whitney 
et al [8], Rutherford et al [9], and The Advanced Composites Design Gu ide 
[10], this test method consists of bolting or bonding stiff rails to a 
composite plate to form a long, narrow test area between the rails. The 
load is introduced at the ends of the rails to displace the rails es- 
sentially parallel to one another. 

A Fourier series plate solution and parametric study of this speci- 
men was conducted by Whitney et al [8]. The study indicated that speci- 
men dimensions can be standardized to produce an acceptable area of 
shear stress for most laminates. The only exceptions are laminates with 
a high Poisson's ratio {= 1.0). Experimental results presented in this 
study showed excellent agreement with lamination theory for in-plane 
shear modulus. A double rail shear test procedure was investigated by 
Sims [11] who found good agreement with both the [±45]^ tension coupon 
and the torsion tube for linear and nonlinear shear stress-strain response. 
Bert [5] pointed out that the rail shear specimen may be readily sized 
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to avoid buckling. 

2.5 The Picture Frame Shear Specimen 

The picture frame specimen shown in Fig. 2c consists of two thin 
square composite laminates bonded to a honeycomb core material for 
stability. A rigid four bar linkage frame is then bolted over the 
sandwiched composite for load introduction. In order to reduce high 
corner stresses these areas usually are cutout and have doublers bonded 
on near the resulting free edge. The load is introduced by a combina- 
tion of tension and/or compression forces in pairs on opposite corners 
Thus, the specimen is subjected to shear loading along planes at 45“ to 
the diagonals. 

A state of pure shear was not observed in a photoelastic study of 
this specimen by Bryan [12]. Dastin et al [13], Hyer and Douglas [14], 
and Hadcock and Whiteside [15] tested boron-epoxy and borsic-aluminum 
picture frame specimens with cutouts and/or doublers at the corners; 
Hadcock and Whiteside [15] and Douglas et al [16] reported that the' 
corner stress concentrations nevertheless remained significantly high in 
both experiments and in the finite element analysis that included 
frame and nonlinear effects. However, Dastin et al [13] reported that 
reliable initial shear modulus and ultimate strength values could be 
obtained from the picture frame specimen for [0/90]^ laminates. 

2.6 The [^45]g Tension Coupon 

The [±45]j tension coupon is a simple test for obtaining linear and 
nonlinear shear stress-stain response in the natural lamina coordinates. 
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As Shown in figure 2d, the specimen is simply a [±45]^ tension coupon, 
usually with tabs. During loading, incremental values of axial force 
and axial and transverse strain are recorded. Various transformation 
laws as employed by Petit [6], Adams and Thomas [17], Rosen [18], and 
Hahn [19] may be used to obtain the piecewise linear modulus of rigid- 
ity, G.|2, in the natural lamina coordinates. 

Petit compared test results of the [±45]^ tension coupon to those 
of the cross beam and found close agreement with the exception of the 
highly nonlinear portion of the shear stress-strain curve. He also 
pointed out that the laminate Poisson's ratio must be exactly 1.0 in 
order to have a state of pure shear in the [±45]^ tensile coupon. The 
measured value of Poisson's ratio for his tests was 0.85. The study 
also reported that specimen delamination may take place due to edge 
effects at high strain values. 

2.7 The Off-axis Unidirectional Tension Coupon 

The off-axis unidirectional tension coupon is a simple specimen 
similar in geometry to the [±45]^ tension coupon shown in figure 2d. 

The laminate is unidirectional with fibers oriented at an angle usually 
between 10 and 75 degrees to the axis of loading. This specimen has 
been used to measure linear and nonlinear shear response in the natural 
lamina coordinates (with the aid of transformtion equations). 

The state of stress in the off-axis unidirectional tension coupon 
has been determined by several methods. Rizzo [20] and Richards et al 
[21] analyzed this specimen by means of anisotropic finite elements; 
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geometry, boundary conditions, and influence of material properties were 
among the parameters studied. Pagano and Halpin [22] described the 
stress state by means of an elasticity solution supplemented with ani- 
sotropic rubber specimens subjected to large deformations. 

Pipes and Cole [23] showed good agreement between the shear re- 
sponses as obtained from 15®, 30®, 45*^, and 60*^ off-axis coupons when 
corrections for Poisson effect are included. Richards et al [21] 
obtained good agreement between experiment and their finite element 
solution using long narrow 45'’ boron-epoxy coupons. Wu and Thomas [24] 
developed grips that allowed in-plane rotations of the anisotropic off- 
axis tension coupon. Chamis and Sinclair [25] suggested an optimum 
fiber orientation of 10° to minimize the effects of the biaxial state of 
stress present in the off-axis coupon. 

2,8 The Torsion Tube 

A desirable specimen for determining shear stress-strain response 
is the torsion tube (Bert [5,26]). In this test a thin-walled circular 
tube is subjected to pure torque resulting in a state of uniform pure 
shear away from the ends. The torsion tube has a very desireable shear 
stress distribution with minimal normal stresses in the test section. 
Torsion tubes have been successfully used to obtain both linear and 
nonlinear shear stress-strain response, however, Rizzo and Vicario [27] 
reported that strength values may be influenced by gripping the ends of 
the tube. The main drawback to this specimen is the difficulty associ- 
ated with fabricating and testing tubular specimens. 
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2.9 Other Shear Test Methods 

Various other test methods have been used to measure in-plane shear 
response of composite materials. See Bert [5], Adams and Thomas [17], 
Hancox [28], Greszczuk [29], Sattar and Kellogg [30], and Berg et al 
[31] for a description of these methods. These other methods generally 
provide qualitative data rather than accurate, quantitative shear 
stress-strain response. 

2.10 Summary 

While this literature review has discussed many test methods for 
the determination of in-plane shear response of composite laminates, it 
is evident that complete analyses including effects of specimen geometry 
and grips, laminate configuration, high temperature test requirements, 
and nonlinear behavior are lacking. This study will attempt to resolve 
some of these deficiencies for flat specimens. 


3. THEORETICAL CONSIDERATIONS 


3.1 Lamination Theory 

This study is limited to linear, elastic, orthotropic laminates in 

plane stress due to mechanical and thermal loads. For such material 

behavior, the laminate moduli of elasticity (E ,E ), Poisson’s ratios 

X y 

(v ,v ) modulus of rigidity (G ) , and the coefficients of thermal 
xy yx xy 

expansion (a ,a ) referenced to the laminate x-y coordinate system 
X y 

(Fig, 4) are the only material properties required for a finite element 
analysis. Experimental results have shown that these elastic and ther- 
mal constants of the laminate can be predicted given the lamina 
properties. The accepted method for doing this is called lamination 
theory [32]. 

The lamina, or individual ply, is the basis upon which the lamin- 
ate is built. It is an orthotropic ply in a state of plane stress and 
is characterized by the constitutive relationship: 


TQ 
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(3.1) 


where the stresses a.j , strains e.j , e^, are referred to 

the 1-2 natural (or principal) lamina coordinate system shown in Fig. 4. 
The terms of the lamina stiffness matrix, [Q], can be expressed in 
terms of the engineering constants as follows: 


<11 


Ei^/(1-v.,,v5J 


12"21 


13 


14 



z 

FIGURE 4. LAMINA 8 LAMINATE COORDINATES 


15 


^22 ~ 

^12 ” ^ 21 ^ 1 l /^^"'' l 2 '' 21 ^ " '' 12 ^ 22 ^^^"^ 12 '' 21 ^ 
^66 " '^12 
^16 *^26 ” 


where all material properties are referred to the lamina principal 
directions. When the principal material directions of any lamina are 
not coincident with the laminate x-y coordinate system, the constitutive 


relations for each lamina must be referred to the laminate x-y coor- 
dinate system. This is accomplished using the transformation matrix 
[T], defined as 


[T] = 


m 


and the Inverse [T]'^ defined as 


m 


2mn 

“2mn 


2 2 

L-mn mn m -n^ 


r 2 

m 


[T] 


-1 


m 


"2mn 

2mn 


2 2 

L -mn m -n . 


( 3 . 3 ) 


where 


m = cose 
n = sine. 


Theta (e), the angle between the lamina natural coordinate system and 


the laminate x-y coordinate system, is shown in Fig. 4. Using the 
transformation matrix to relate the stresses and strains in the natural 
coordinate system to those In the x-y coordinate system leads to the 
relationship 

‘"x'k ' 

Thus, for the lamina, the stresses, >"elated to the 

strains, {e }. , in the x-y coordinate system by the transformed lamina 

X K 

stiffness matrix, [Q]|^. defined as 

[Q]|^ = [T]-’[qim,,. (3.5) 

For symmetric laminates subjected to in-plane loading only, strains in 
each lamina are equal to the laminate midplane strain, {e°}: 

= {e°}. (3.6) 

The stress resultants per unit length, {N>, are defined in terms of the 
stresses by the relationship 



where 2H is the total laminate thickness. Combining equations (3.4), 
(3.6), and (3.7) and rewriting the integral in terms of a summation 
over all n plies, the force resultants are 
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From equation (3.8), the extensional stiffness matrix Is defined 


so that the force resultant relationship can be written 

{N} - CA]{£®}. 

The average laminate stresses, {a}, are defined 

1 


(3.9) 


(3.10) 


(ct) 


2H 


(N). 


(3.11) 


Combining equations (3.10) and (3.11) and considering appropriate non- 
zero {a} terms, the effective laminate engineering constants are: 

E. 


^ 11 ^ 22'^12 


“2M22 

A A -A^ 

nr22 ”12 




xy 


yx 


11 

^12 

^22 

Ai2 


(3.12) 


P . ^66 

^xy " 2H • 

The laminate thermal coefficients of expansion, {a}, are derived 
in a similar manner. When the lamina thermal strains, 
included the lamina constitutive equation (3.4) becomes 




(3.13) 
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and with zero force resultants, equation (3.8) becomes 

0= ” (3.14) 

1 

Solving (3.14) for the laminate thermal strains {e^} results in 

= [A]-^ " 
k” 1 

and for uniform temperature change aT, the lamina thermal strains are 

{cjlk = 

Now defining the laminate thermal coefficients {a} as 

(3.17) 

and combining (3.15), (3.16), and (3.17) gives the relationship 

{5} = [A]”’’ Z CQ],^{ax}k(h,^-h^_.,). * (3.18) 

k“ 1 

3.2 The Finite Element Formulation 

The finite element method is a versatile technique for obtaining 
approximate solutions for structures subjected to loads. Text book 
accounts can be found in Desai and Abel [33] and Zienkiewicz [34], 
among others. This section describes the theory of the two-dimensional 
finite element program used in this study. The analysis considers 
orthotropic materials in plane stress due to mechanical and thermal 
loading. 

The program uses both triangular and quadrilateral elements with 
the quadrilateral elements being composed of four triangles. All 
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triangles are constant strain elements. Convergence requirements for 
the finite element solution are satisfied by using a line u' displacement 
model for the in-plane displacements U and v in the x and y directions, 
respectively; 


u = Bi+82X+p3y 


V = g 4 + 65 X+Bgy, 


(3.1SI 


where the generalized coordinate coefficients, , are functions of the 
initial location of generic nodes 1, 2 , and 3 and are given by 


where 



^2^3'^3^2 

^3 "^2 

X3-X2 


ys-yi 

X3-X1 

_^1'5'2"Vi 

y 2 -y-| 

X2-X1 

^2'''3‘^3^2 

^3-^2 

X 3 -X 2 

^1^3-^! 

^3-^1 

X3-X1 

^ 1 ^ 2 “Vl 

yg-yi 

X2-X1 




The general strain-displacement relationships are 


(3.20) 


p B 

X ax 
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3V 

9y 


= 9U . 3V 

y es ■ ' ♦f' 

'xy 9y 9x 


(3.21) 


and combining equations (3.20) and (3.21) leads to the specific strain- 
displacement relationship 

{£} = [B]{u> (3.22) 


where the strain-displacement matrix, [B], and the transposed nodal 
displacement matrix, {u}^, are defined 


[ 6 ] = k 

>2-^3 

^3-^1 

0 

CVJ 

>r 

0 

0 


0 

0 

0 X 3 -X 2 

X 1 -X 3 

x^-Xi 

(3.23) 

and 

X 3 -X 2 

X 1 -X 3 
[u, Uj 

X 2 -X 1 y2-y3 

“3 '^l ^^2 

^3-^1 

V 3 ], 

^ 1-^2 

(3.24) 


where A is the area of the element. 

The thermoelastic constitutive relationship for an orthotropic 
laminate is 

{a} = [a]({e}-{e^), (3.25) 

where the strain that contributes to the laminate stress is the dif- 
ference between the total midplane strain, {el, and the strain caused 
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V 


only by initial (free thermal) strains, (eq. 3.15). The terms of 
the orthotropic laminate stiffness matrix, [a], referenced to the x-y 
coordinate system, are defined 

*11 “ V<'-V.vx> 

*22 ° 

*12 " ' %yV<'"'’xy'’yx’ 

*66°*^xy 
*16 " *26 ' “• 

The total potential energy, n, of a system is the sum of the strain 
energy, U, and the potential energy of the external loads, W: 

n = U + W. (3.27) 

Considering equations (3.22-3.25), U and W are, respectively 

U = I fff {e}''’[a]{e}dv - /// {e}^[a]{E^}dv 
V V 

+ i/// {eb'^[a3{eSdv (3.28) 

^ V 

W = - {u)’^{Q} 

where V is the elemental volume and the transposed mechanical nodal load 
vector, {Q}^, is defined as 

= [q^ d2 i3 P-| p2 Pa^ 


(3.29) 


22 


The variables and p^. are mechanical loads applied in the x and y 
directions, respectively, at the i^*^ node of the element. Minimization 
of the total potential energy of an element (eq. 3.27) leads to the 
equilibrium equations 

[k]{u) = {Q} + {qb (3.30) 

where the relationships for the elemental stiffness matrix, [k], and the 
initial load vector, {Q^}, are 

[k] ^ fff [B]’^[a][B]dv (3,31) 

V 

= /// [B]"'’[a]{0''}dv (3.32) 

V 

The total stiffness matrix [K] and the total load vector {R} for the 
entire structure is formed by combining the elemental stiffness matrices 
and the total load vector ({Q}+{Q^}) for all the elements resulting in 
the system of linear equations 

[K]{X> = {R}. (3.33) 

For N number of nodes the total stiffness matrix will be 2N x 2N and the 
total load vector and node displacement vector {X} will be 2N. These 
equations relate the stiffnesses and nodal displacements of the entire 
structure to the applied loads and hence this formulation is called the 
direct stiffness method. Geometric boundary conditions are then imposed 
and a solution is obtained for the node displacement vector, {XI. 


4. THE SLOTTED SHEAR SPECIMEN 


4.1 General 

The slotted shear specimen was analyzed using the basic finite 
element representation shown in Fig. 5. The specimen was always loaded 
to simulate uniform axial and zero lateral displacement at the ends of 
the specimen. These boundary conditions are idealizations of machine 
grips loading a flat coupon specimen in tension. The effects of lam- 
inate properties were studied for three graphite-epoxy laminates. The 
[0], [90], and [±45]^ laminates were chosen because they offer sub- 
stantial variations in orthotropic properties. The [0] laminate ex- 
hibits a low value; the [90] laminate exhibits a high E^/E^ 

value; and the [±45]^ laminate exhibits a high G^y/Ey value. 

Elementary rigid body statics assumes the shear stress distributioini 
in the test section between the ends of the slots to be pure, uniform, 
and maximum and independent of the laminate properties. Actual stress 
contour plots for the laminates studied are shown in Figs. 6-8. All 
stresses have been normalized with respect to the average axial stress 
away from the test section. Since a limited parametric study on the 
slot spacing, S (ranging from 0.2W to W), and slot depth, D (ranging 
from 0.25W to 0.75W), showed no significant changes in the stress 
distributions, the slot spacing, slot depth, and slot width (T) were 
standardized to S = 0.6W, D = 0.5W, and T = 0.0625W for all subsequent 
studies. As indicated in the figures, the finite element results did 
exhibit the expected symmetries. 
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FIGURE 5. FINITE ELEMENT REPRESENTATION 
OF SLOTTED SHEAR SPECIMEN 
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4.2 [0] Laminate 

The normalized stress contours for; the [0] slotted shear specimen 

are shown in Fig. 6. The a /a and a /a contour plots clearly indicate 

X y 

stress concentrations at the slot ends (a /o = 0.4 and a /a = 8.0) 

X y 

with the Oy/a stress remaining at 4.0 along the centerline of the test 
section. The a /c stress is insignificant (o /a = -0.065) at the 

^ X 

center of the specimen. The normalized shear stress remains essentially 
constant throughout the test section at -0.6; the maximum magnitude is 
-0.8 near the slot ends. Although this stress distribution is ac- 
ceptable in principle, this region may be too narrow for practical 
strain measurements. 

Since a region of uniform shear stress occurs in the test section, 
the [0] slotted specimen may be used for qualitative shear measurements. 
However, stress concentrations at the slot ends and a rather small area 
of uniform shear stress make this specimen unacceptable for quantitative 
shear measurements. 

4.3 [90] and [±45]^ Laminates 

The contour plots for [90] and [±45]^ graphite-epoxy laminates are 

shown in Figs. 7 and 8. The stress distributions for these laminates 

are clearly unacceptable for a shear specimen. As well as having 

a^/o and 0 y/a stress concentrations at or near the slot ends, the 

T^y/a contour plots indicate that the state of shear is not maximum or 

uniform (for the [90] laminate, t Jo ranges from -2.0 to 0.25 and for 

xy 

the [±45] laminate, t /o ranges from-l.C to -0.4). 

^ xy 
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FIGURE 6. NORMALIZED STRESS CONTOURS FOR [O] 
SLOTTED SHEAR SPECIMEN ,S=0.6W,D=05W,T=0.0625W 
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FIGURE 7 NORMALIZED STRESS CONTOURS FOR [ 90 ] 
SLOTTED SHEAR SPECIMEN , S=0.6W.D=0.5W,T=0.0625W 
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FIGURE 8. NORMALIZED STRESS CONTOURS FOR [±45 
SLOTTED SHEAR SPECIMEN, S=0.6W,D=05W,T=0D625' 
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4.4 Summary 

For the [0], [90], and [±45]^ laminates, the slotted shear specimen 
exhibits 0 , o , and stress concentrations in the vicinity of the 
slot ends as expected. As a result of such stress concentrations, high 
normal stresses were often found in the test section. For the [90] and 
[±45]^ laminates, is very low in magnitude (for the [90] laminate 

T /o = 0.25 and for the [±45] laminate t /o = -0.4) and nonuniform in 

xy 

the test section. Obviously this slotted specimen does not provide 
acceptable shear data for these laminates. The normalized shear stress 
in the test section of the [0] laminate is somewhat uniform (with the 
exception of a small area in the Vicinity of each slot end), however, 
this area of high shear stress is rather narrow for strain measurements. 
Also, the normalized 0 ^ is significant in the test section for the [0] 
laminate. In conclusion, the stress distributions* of the slotted shear 
specimens analyzed were significantly influenced by laminate properties 
and none of the laminates studied were found to be acceptable for 
providing reliable quantitative shear data. 


5. THE CROSS BEAM SHEAR SPECIMEN 


5.1 General 

As mentioned previously, the cross beam shear specimen shown in 
Fig. 2a has been used to obtain shear stress-strain response for both 
[±45]g and [0/90]^ laminates. When analyzing the results from such a 
test it is assumed that the test section is subjected to uniform biaxial 
normal stresses of equal magnitude and opposite sign resulting from the 
application of positive and negative bending moments to the cross beam. 
For this state of biaxial stress, stresses on 45° planes are pure shear 
and equal to the magnitude of the applied average normal stresses. The 
purpose of this analysis is to study the influence of corner effects, 
laminate properties, and, to a limited extent, core materials on the 
stress distribution in the test section. 

The elastic analysis of the cross beam shear specimen involved the 
following assumptions. The composite material was assumed to be in a 
state of plane stress. This neglected any transverse stress components 
associated with the application of loads near the ends of the beam (Fig. 
2a) and also interlaminar effects from the core in the test section. 
Loading and geometric symmetry allowed the analysis to be applied to 
only one quarter of the central portion of the specimen. These symmetry 
conditions were simulated by imposing zero x displacement along the y 
axis and zero y displacement along the x axis shown in Fig. 9. 

Cross beam specimens usually are designed with two types of honey- 
comb core. Stiff, heavyweight core material is sandwiched between the 
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FIGURE 9. FINITE ELEMENT REPRESENTATION 
OF CROSS BEAM SHEAR SPECIMEN 
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composite material away from the central portion of the specimen (x and 
y > W) to preclude local failure at the points of load application, 
while a more flexible, relatively light weight core material is used in 
the central portion (x and y < W). Because of the unknown net stiffen- 
ing effects of core, displacements rather than forces were imposed as 
boundary conditions at x = W md y = W to represent the loading intro- 
duced by the bending moments. In order to study core stiffness effects 
away from the test section, two sets of lateral displacement boundary 
conditions were considered at x = W and y = W. "Rigid" core for x and 
y greater than W was simulated by fixing lateral displacements to zero 
across the width of the beam. "Flexible" core was simulated by allowing 
free lateral displacements across the width of the beam. The stiffness 
of the core material in the test section and all other core influences 
were always neglected. 

Contour plots of the stresses on a 45'’ plane (x',y'), normalized 
with respect to the magnitude of the applied average normal stress, o, 

are shown in Figs. 10-14. These plots reflect the symmetric loading and 
geometry. 

5.2 The [±45]^ Laminate 

A cross beam with fibers at +45“ is used to determine shear stress- 
strain response of [0/90]^ laminates. Such response is equivalent to 
the shear response of a lamina referenced to its natural coordinates. 
Figures 10-12 compare the rigid core results to the flexible core 
results for normalized o', c' and t' stress distributions. 

A y xy 
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For both cases the shear stress is maximum and pure over a finite 

region in the center of the specimen. This region is large enough for 

a strain gage to be used for accurate strain measurements. Although 

these are excellent features of the cross beam, the magnitude of 

t' /a at the center is significantly different from 1.0 because of the 
xy 

corner effects. Specifically, values of t'Jo were found to be 1.20 and 
1.13 for the rigid and flexible core boundary conditions, respectively. 

Corner stresses for the rigid core boundary conditions were lower 
than those of the flexible core boundary conditions, especially for the 
Q^/a stress. With one exception (Fig. 11a), all normalized stress com- 
ponents are greater than 0.6 in the vicinicy of the corner. This is 
significant and may contribute to premature failure in a mode different 
from pure shear. 

In general the stress distributions in the [±45]^ cross beam are 
significantly more complex than those of the "ideal" cross beam. The 
magnitude of the shear stress in the test section is higher than 1.0 due 
to corner effects and stress concentrations at the corners undoubtedly 
contribute to any failures. However, this analysis indicates that the 
stress distribution in the center section is uniform pure shear. Thus, 
the elastic shear modulus may be determined knowing the magnitude of the 
shear stress in this region. 

5.3 The [0/90] j Laminate 

A cross beam with fibers at O'” and 90° Is used to determine shear 
stress-strain response of a [±45]^ laminate. Finite element results for 
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FIGURE 13. NORMALIZED S a-' CONTOURS 
FOR [0/90]g CROSS BEAM 
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FIGURE 14. NORMALIZED CONTOURS 
FOR [0/90]g CROSS BEAM 
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both the rigid and flexible core boundary conditions were similar and 
thus specific results for both core types are shown only for the 
normalized stresses. The results are not affected appreciably by 
the lateral displacement boundary condition because the high lateral 
stiffness of the specimen results in very small lateral displacements. 

As in the [±d5]^ laminate, the [0/90]^ laminate produces a maximum, 
pure shear region in the center of the specimen. This shear region has 
more characteristics of the “idear shear specimen. An area of maximum, 
uniform shear stress is present over nearly the entire test section and 
the maximum is very close to 1.00 (exactly 1.01). In the vicinity 

of the corner of the specimen all the normalized components of stress 
remain above 0.6, however. 

The finite element analysis indicates that the [0/90]^ laminate 
should provide reliable stress-strain response to a high degree of 
accuracy in the linear portion of the shear stress-strain curve. Al- 
though failures may initiate due to corner stresses, the core effects 
are minimal with the [0/90]^ laminate. 

5.4 Summary 

For both the [±45]^ and the [0/90]^ cross beam shear specimens, the 
maximum value of occurs at the center of the specimen and is suit- 
able for strain measurements. However, the magnitude of x* /o at the 
center is dependent on the laminate properties, and for the [±45] 
laminate, x^a is significantly higher than 1.0. For laminates that 
are relatively stiff across the width of the beam compared to their 
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Shear stiffness, the honeycomb core effects were shown to be minimal 
with this plane stress analysis. 

Finite element results indicate the cross beam may be used with a 
[0/90]^ laminate for linear response since its stress distribution is 
nearly "ideal" (with the exception of high corner stresses). The state 
of shear stress at the center of the laminate is influenced by 

the core stiffness and is significantly greater than 1.0 because of a 
low shear stress away from the center. Thus, the applied normal stress 
in the [±45]^ laminate is not simply related to the shear stress at the 
center of the specimen. However, initial shear modulus values may be 
determined if this relationship is known from a numerical study such as 
this one. 


6. THE lOSiPESCU SPECIMEN 


6.1 General 

The specimen analyzed in this study is a modified losipescu specimen 
[3] in that the specimen is a flat laminate notched on two sides with 
the load introduced through fixtures bonded to the specimen. losipe^cu's 
original specimen was a bar notched on four sides and loads were applied 
directly to the specimen. This modified specimen was studied with both 
isotropic and orthotropic materials. In addition, several minor varia- 
tions in the specimen configuration were studied as were the influence 
of elastic fixtures and thermal loads. 

6.2 Rigid Fixture Results 

6.2.1 General 

Several laminates were studied assuming that the fixtures are rigid 
(Fig. 2b). The finite element mesh used to represent these specimens is 
shown in Fig. 15. Because a finite element analysis of the entire 
specimen indicated quarter symmetry of stresses and displacements, only 
one quarter of the specimen was modeled in subsequent studies. This was 
accomplished by applying uniform u and zero v displacements along the 
line y = L/2, and constraining the v displacements along x = 0 and u 
displacements along y = 0 as indicated in Fig. 15. The dimensions of 
the specimen were standardized to: notch depth, A = 0.225W; notch angle, 

0 = 90°; and specimen length, L = W, where W is the width. Normalized 
stress contour plots for the rigid fixture assumption are shown in Figs. 
16-21 and results are tabulated in Table 1. All stresses are normalized 
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FIGURE 15. FINITE ELEMENT MODEL FOR 
lOSIPESCU SPECIMEN WITH RIGID FIXTURES 
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with respect to t, the average shear stress along the line OA between 
the notches (Fig. 15). 

6.2.2 Isotropic and Quasi-isotropic Materials 

Normalized stress contours for an isotropic material (steel) and a 
quasi-isotropic material (C0/9O/±45]^ graphite-polyimide) are shown in 
Figs. 16 and 17, respectively. The two stress distributions are nearly 
identical exhibiting essentially uniform pure shear in the test section. 
The normalized shear stress along the centerline OA ranges from 0.97 at 
the center 0 to 1.12 in the vicinity of the notch tip A (Table 1). As 
expected, stress concentrations are present at the corner C. These 
linear elastic results indicate that this specimen produces a very 
satisfactory stress distribution for determining the linear shear re- 
sponse of isotropic and quasi -isotropic materials. They are in agree- 
ment with the photoelastic results obtained by losipescu. It is also 
evident from these results that failure will initiate at the notch tip A 
or the corner C due to the stress concentrations in these regions. 

6.2.3 The [0] Laminate 

Normalized stress contours for a [0] laminate are shown in Fig. 18. 
The shear stress is pure, uniform, and maximum at the center of the 
specimen (point 0), but is not uniform along the centerline OA, ranging 
from 1.10 at the center 0 to 0.48 at the notch tip. Stress concentra- 
tions are again present at the corner C. These linear elastic results 
indicate that the stress distribution at the center of the [0] specimen 
is suitable for linear shear response, however, the magnitude of this 
normalized shear stress is greater than 1.0 due to the nonuniform shear 
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FIGURE 16. NORMALIZED STRESS CONTOURS 
FOR lOSIPESCU SPECIMEN WITH ISOTROPIC 
MATERIAL AND RIGID FIXTURES 
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FIGURE 17. NORMALIZED STRESS CONTOURS FOR 
[0/90/t45]s GRAPHITE- POLYIMIDE lOSIPESCU 
SPECIMEN WITH RIGID RXTURES 
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stress distribution along the centerline OA. Since all stresses were 
rather low in magnitude at the notch tip, failure is expected at the 
corner C . 

6.2.4 The [90] Laminate 

Normalized stress contours for a [90] laminate are shown in Fig. 19. 
Although a pure, uniform shear stress at the center 0 is present, the 
distribution is very nonuniform along the centerline OA with norma- 
lized shear stress values ranging from 0.84 at the center to 1.94 at the 
notch tip A. Because this laminate exhibits a nonuniform shear stress 
distribution along OA, a stress factor relating the shear stress at the 
center of the specimen to the applied load (0.84) must be used in order 
to obtain accurate initial shear response. In addition, stress concen- 
trations at the notch tip A will contribute to any failures. 

6.2.5 The [0/90]^ Laminate 

Normalized stress contours for a [0/90]^ laminate are shown in Fig. 
20. The normalized shear stress along the centerline OA ranges from a 
maximum of 1.09 near the notch tip to a value of 0.96 at the center. As 
expected, stress concentrations are again present at the corner C. Since 
this laminate produces a very uniform region of pure shear in the test 
section, it will provide accurate linear shear response. These results 
also show that failure may initiate at the notch tip A or at the corner 
C due to the stress concentrations in these regions. 

6.2.6 The [±45]^ Laminate 

s 

Normalized stress contours for a [±45]^ laminate are shown in Fig. 




FIGURE 18. NORMALIZED STRESS CONTOURS 
FOR [0] GRAPHITE - POLYIMIDE lOSIPESCU 
SPECIMEN WITH RIGID FIXTURES 
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FIGURE 19. NORMALIZED STRESS CONTOURS 
FOR [90] 6RAPHITE-P0LYIMIDE lOSIPESCU 
SPECIMEN WITH RIGID FIXTURES 
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FIGURE 20. NORMALIZED STRESS CONTOURS 
FOR [0/90]s GRAPHITE- POLYIMIDE lOSiPESCU 
SPECIMEN WITH RIGID FIXTURES 



0 . 




FIGURE 21. NORMALIZED STRESS CONTOURS 
FOR [t45] GRAPHITE- POLYIMIDE lOSIPESCU 
SPECIMEN WITH RIGID FIXTURES 
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TABLE 1. Normalized Stresses for Graphite-polyitnide and Steel 
losipescu Specimen with Rigid Fixtures. 



oi 

a1 

cl 

Steel 

ox/t 

e.z/xio-'^ 

0.42 

0,89 



ay/V 

4.18x10'5 

0.35 



0.97 

1.05 

0.94 

[0/90/±45]s 

Gr/Pi 

ax/t 

6.28x10-4 

0.42 

0.82 

ay/7 

4.08x10-5 

0.35 

4.61 

•rxy/7 

0.97 

1.05 

0.88 

[0] 

Sr/Pi 

Ox/7 

3.41x10-3 

0.25 

0.57 

Oy/x 

2.12x10-4 

0.18 

3.26 

^xy/^ 

1 .10 

0.48 

0.43 

[90] 
Sr/ Pi 

o^xA 

3.34x10-4 

0.56 

8.46x10-2 

ay/T 


0.62 

6.91 

Txy/x 

0.84 

1.95 

0.38 

[0/90]s 

Gr/Pi 

Ox/t 

2.06x10-3 

0.43 

0.13 

ay/T 

-7.92x10-5 

0.34 

6.26 

Xxy/7 

0.96 

1.02 

0.40 

[±45]s 

Sr/Pi 

ax/t 

2.38x10-4 

0.41 

3.32 

ilHEHl 

7.18x10-4 

0.33 

6.38 

_ 

Txy/x 

1.09 

1.06 

3.17 


^See Figure 15. 
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21. The shear stress distribution in the vicinity of the center of the 
specimen is essentially pure and uniform. The normalized shear stress 
along the centerline OA ranges from a maximum of 1.10 near the notch tip 
to a minimum of 0.89 between 0 and A. The value at the center of the 
specimen is 1.09. All stress components are unusually high at the 
corner C. These linear elastic results indicate that this specimen 
produces a very satisfactory stress distribution for determining the 
linear shear response of a [±45]^ laminate. It is also evident that 
failure will initiate at the corner C due to the stress concentrations 
in this region. 

6.3 Elastic Fixture Results 

6.3.1 Mechanical Loading 

The finite element representation for the specimen with elastic 
fixtures is shown in Fig. 22. This finite element representation is a 
two-dimensional elastic model of a specimen with four steel fixtures 
(each with a thickness T^ = 0.25W) bonded to the test material. Lamina- 
tion theory was used to predict the elastic properties of the fixture - 
composite laminate. All dimensions were the same as those used in sec- 
tion 6.2 and in addition; A = 0.25W; fixture dimensions, F-| = 0.65W, 

F^ = 0.25W, Fj = O.IW; fixture angle, = 45°; and specimen thickness, 
T^ = 0.0384W ( 10 plies of graph! te-polyimide, 0.00768"/ply) . The 
stresses were normalized with respect to the average shear stress deter- 
mined from the given applied load and the area of the specimen between 
the notch tips A and B. 
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FIGURE 22. FINITE ELEMENT MODEL FOR 
lOSIPESCU SPECIMEN WITH ELASTIC FIXTURES 
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6. 3. 1.1 Isotropic Material 

Normalized stress contours for an isotropic specimen (steel) are 
shown in Fig. 23. A pure uniform shear stress distribution exists along 
the centerline AB with normalized values ranging from 0.89 at the notch 
tips A and B to a maximum value of 0.93 at the center of the specimen. 

It is obvious that these results, which were obtained from a rather 
coarse finite element mesh, do not satisfy overall force equilibrium 
exactly since the total shear force transmitted across the test section 
is not equal to the applied load. 

These results can nevertheless be compared to the rigid fixture 
solution (Fig. 16) in a general manner. With the exception of the lack 
of quarter symmetry near the elastic fixtures, the shape of the stress 
contours are very similar. The largest difference occurs near the corners 
C-F of Fig. 23 where the deformation of the elastic fixtures has the 
greatest effect on the specimen. Along the centerline AB (with rigid 
or elastic fixtures) the effects of bending in the fixtures are clearly 
diminished because of the symmetries that must exist there. This feature 
of the stress distributions is evident when the shapes of the 0.6 and 
0.9 contour lines in Fig. 23c are compared. Therefore, although failure 
may initiate at the points A-F, this specimen produces a satisfactory 
stress distribution for determining linear shear response of isotropic 
materials. 

6. 3. 1.2 The [0] Laminate 

The normalized stress contours for a [0] laminate are shown in Fig. 

24. The shear stress distribution along the centerline AB is essentially 




FIGURE 23. NORMALIZED STRESS CONTOURS FOR lOSIPESCU 
SPECIMEN WITH ISOTROPIC MATERIAL AND ELASTIC FIXTURES 
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pure, uniform, and maximum at the center of the specimen. The normalized 
values range from 0.72 at the notch tips A and B to 1.03 at the center 
of the specimen. Similar to the isotropic specimen, the fixture effects 
are minimal as a comparison of Figs. 18 and 24 shows. It is also 
evident from these results that failure is expected to initiate at the 
notch tips A or B or at the corners C, 0, E, or F due to the stress 
concentrations in these regions. 

6.3.2 Thermal Loading 

In this section the effects of applying a temperature change of 

-600°F to a complete stress free specimen (laminate bonded to fixtures) 

are studied. This loading simulates the elastic behavior resulting from 

an elevated temperature bond of the fixtures to the specimen. The [0] 

laminate was used for this analysis with the mesh shown in Fig. 22 (F = 

0). Contour plots for a , a , and t normalized with respect to their 

A y xy 

corresponding ultimate stress values, (f^^^ult’ ^^y^ult’ ^"^xy^ult 
are shown in Fig. 25. 

The most significant feature of the results is that the thermal 
stresses are negligible in the test section for all stress components. 
Normalized a and t stresses are both very high at the corners with 

y ^y 

maximum values of -2.49 and 1.72, respectively. As expected, the stress 
in the fiber direction was quite low throughout the specimen with a 
maximum normalized value of -0.29 at the fixture-specimen interfaces CD 
and EF. If the specimen was mechanically loaded and the behavior was 
elastic, a superposition of the stresses depicted in Figs. 24 and 25 


57 


% 



FIGURE 24. NORMALIZED STRESS CONTOURS FOR [0] GRAPHITE - 
POLYIMIDE lOSIPESCU SPECIMEN WITH ELASTIC FIXTURES 


AT=-600®F 



FIGUF^ 25. NORMALIZED STRESS CONTOURS FOR THERMALLY LOADED [C3 
GRAPHITE- POLYIMIDE lOSIPESCU SPECIMEN WITH ELASTIC FIXTURES 
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would give the combined stress distribution. Although the corner 
stresses at D and E would cancel one another for combined mechanical and 
thermal loading, the stresses at corners C and F would be additive. For 
this reason fixtures should be bonded to specimens at the desired test 
temperature in order to minimize the stress concentration at corners C 
and F. 

6.4 Other Results 

6.4.1 Doubler Effects 

In addition to the elastic and rigid fixture cases described pre- 
viously, the effect of doubling the stiffness (doubling the thickness in 
practice) of the test specimen in the vicinity of the corners was studied 
The results showed that other things being equal (dimensions, loads, 
etc.), doubling the stiffness in the vicinity of the fixture/specimen 
interface reduced the stress concentrations in this region by approxi- 
mately 50 percent (Table 2) while the stress distribution in the center 
of the specimen remained virtually unaltered. These results are signi- 
ficant because they indicate that the high corner stresses can be re- 
duced to reasonable values through the use of doublers without signi- 
ficantly effecting the stress distribution in the test section. 

6.4.2 Aspect Ratio Effects 

The influence of specimen aspect ratio, R = L/W (Fig. 22), was 
investigated for R = 1.0 and R = 2.5. The stresses in the test section 
of the specimen were not appreciably effected for the elastic and rigid 
fixture cases studied, however, the corner stresses were significantly 


TABLE 2. Effects of Doublers and Aspect Ratios for the losipescu Specimen. 
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TABLE 3. Graphite-epoxy and Graphite-polyimide Normalized Stresses 
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higher for the specimens with R = 2-5 (Table 2). These higher stresses 
were the result of the higher bending moments introduced from the fix- 
tures into the specimen for R = 2.5. The effect of different materials 
(graphite-epoxy and graphite-polyimide) was also studied for R = 2.5 
(Table 3) using a coarse mesh with rigid fixtures. The stress values 
for the graphite-epoxy and graphite-polyimide were virtually identical 
as indicated in the table. With the exception of the [±45]^ laminate, 
the shear stresses for R = 2.5 for these two materials also compare to 
the shear stresses at the center of the specimen for R = 1.0 (Table 2). 

6.4.3 Rounded Notch Effects 

The effects of rounding the notch tips were studied for both [0] 
[90] laminates. This was done by adding a "fillet", or round with a 
radius of 0.0325W to the notch tip. The normalized stress distributions 
along the line y = 0 in the vicinity of the notch tip are shown in Figs. 
26 and 27 for the [0] and [90] laminates, respectively. For the [0] 
laminate all components of stress are considerably lower near the free 
edge with the rounded notch. However, the round causes the shear dis- 
tribution to become more nonuniform, and as a result the normalized 
shear stress at the center of the specimen is higher than that of the V- 
notched specimen (Table 4). As indicated in Fig. 27, rounding the notch 
for a [90] laminate results in a and x approaching zero near the 
free edge. However, o remains very high. Since the shear stress 
distribution for the rounded notch is more uniform, the corresponding 
normalized shear stress at the center of the specimen is more nearly 

equal to 1,0. Thus while a and t are always zero at the free edge of 

X xy 


O^Z .036 .030 .024 .018 .012 .006 NOTCH 

TIP 

s/ w 

FIGURE 26. ROUNDED NOTCH EFFECTS FOR 
THE [0] lOSIPESCU SPECIMEN ALONG Y=0 









042 .036 .030 .024 ,018 012 .008 NOTCH 

TIP 

S/ W 

FIGURE 27. ROUNDED NOTCH EFFECTS FOR 
THE [90] lOSIPESCU SPECIMEN ALONG Y=0 
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a rounded notch, the uniformity of the shear stress distribution is a 
function of the modulus ratio E/Ey- 

6.5 Summary 

While the losipescu specimen produces a region of uniform shear 
stress at the center of the specimen for isotropic and orthotropic 
materials, the uniformity of the shear stress between the notch tips is 
highly dependent upon the elastic properties of the laminate. This 
feature was most evident for the unidirectional laminates. Also, the 
complex state of stress at the notch tips and corners undoubtedly 
contributes to failure. This specimen can accurately provide initial 
shear modulus values for isotropic and orthotropic materials. Bonding 
suitable doublers to the corners and rounding the notch tips can reduce 
the stress concentrations in critical areas. Thermal stresses are not 
significant in the center of the specimen. 


7. THE RAIL SHEAR SPECIMEN 


7.1 General 

As mentioned previously, the rail shear specimen has been used to 
obtain linear and nonlinear shear stress-strain response for composite 
laminates at both room temperature and elevated temperatures. For 
experimental work, it has been assumed that the specimen is subjected to 
pure, uniform shear stress with the magnitude of the shear stress 
directly related to the applied load through a simple P/A relationship. 
The finite element results presented here show to what extent these 
assumptions are valid. The effects of material properties, thermal 
loads, and rail configurations for both rigid and elastic rails are 
presented to more fully understand the elastic behavior of this specimen 
and to more accurately determine the relationship between the applied 
load and the stress distribution in the test section. 

7.2 Rigid Rail Solutions 

7.2.1 Fourier Series Representation 

A Fourier series solution was presented [8] as the superposition of 
two problems with the boundary conditions (Fig. 29} for each problem 
being : 

problem 1. u(x,0) = v(x,0) = v(x,b) = 0 
Nj^(0,y) = N^(a,y) = 0 
u(x,b) = S 

Nx^(0.y) . Nx^(a.y) - f 
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problem 2. u(x,0) = u(x,b) = v(x,0) = v(x,b) = 0 

(0,y) = N^(a,y) = 0 (7.2) 

(0,y) = Nj^y(a,y) = -Agg f , 

where N. and are defined from equations (3.7) and (3.9). These 

1 DO 

boundary conditions assume that the actual specimen is perfectly bonded 
to rigid rails and that uniform rail spacing is maintained during 
loading. The solution to problem 1 leads to a state of pure shear with 
the displacements 

u = ; V = 0. (7.3) 

The solution to the second problem was assumed to be of the form 


u 


m=l n=0 

(o<x<a; o<y<b) 


00 00 


_ r, (2m-l)irx ___ 2mr.y 

2 ^ cos e — cos — ff- 

m=l n=0 


(o<x<a; o<y<b) 


(7.4) 


The superposition of the two solutions result in an infinite m+n system 
of linear equations with unknown Fourier coefficients. Adequate conver- 
gence of the infinite series was obtained by truncating the system of 
equations to 24 terms. The Fourier coefficients were computed and 
stresses were determined as a function of position. 
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7.2.2 Comparison of Fourier Series and Finite Element Results 

Fourier series solutions are compared to finite element results for 
[0/±45]g and [±45]^ Thorne! 50 graphite-epoxy in Fig. 28. The finite 
element representations used are shown in Fig. 29 along with boundary 
conditions for the problem. Results were obtained for a specimen 
aspect ratio a/b = 12. Fig. 28 shows the distribution of the shear 
stress, normalized with respect to the average shear stress, along the 
centerline of the specimen y = b/2. 

As indicated in Fig. 28a, the two solutions exhibited very similar 
results for the [0/±45]g laminate. The normalized shear stress is uni- 
form over almost 80 percent of the specimen at a value of 1.05. This 
value is necessarily greater than 1.0 to account for the shear stress 
approaching zero at the free edge. 

The results for the [±45]^ laminate shown in Fig. 28b do not ex- 
hibit the same high degree of correlation. The finite element results 
are very similar to those of the [0/±45]g laminate with the normalized 
shear stress uniform over almost 90 percent of the specimen at a value 
of 1.06. The Fourier series results exhibit a spike with a maximum 
value of 1.5 a small distance in from the free edge and then exhibit a 
continuously decreasing stress magnitude toward the center of the 
specimen. Because of the difference in these results, additional finite 
element results were obtained using the model shown in Fig. 29b. The 
results obtained with the refined mesh were so similar to those of the 
basic mesh that they cannot be distinguished in Fig. 28b. The finite 
element results are believed to be more accurate than the Fourier series 
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FIGURE 28. C0MR1\RIS0N OF RIGID RAIL SOLUTIONS 
FOR THORNEL 50 GRAPHITE- EPOXY AT Y=b/2 







FIGURE 29. FINITE ELEMENT MODELS FOR 
RAIL SHEAR SPECIMEN WITH RIGID RAILS 
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results^ for this [±45]^ laminate due to the proven capabilities of the 
method, the consistency with the [0/±45]^ laminate, and the more phy- 
sically appealing stress distribution. 

The finite element results indicated that both laminates exhibited 
high normal stresses in the vicinity of corners, but these normal 
stresses were negligible over 80 percent of the specimen. Thus, these 
results indicate that this specimen produces an essentially pure, uni- 
form distribution of shear stress in the test section with the stress 
magnitude equal 1.05 and 1.06 times the average shear stress t for 
[0/±45]g and [±45]^ laminates, respectively. 

7.3 Elastic Rails 

7.3.1 Effects of Offset Loading and Tapered Rails 

The effects of nontapered and tapered rail configurations loaded 
along the centerline, and nontapered rails with offset loading for a 
[0] graphite-polyimide laminate were studied. The finite element mesh 
for the nontapered rail configuration with axial loading is shown in 
Fig. 30. The dimensions of the specimen were standardized to: specimen 

length, L = lOW; fixture dimensions, = 5.83W, F 2 = 3.33W, F^ = 6.67W, 
F^ = 1.23W; rail thickness, T^ = 3.33W; and composite thickness, T^ = 
0.256W. Appropriate steel, steel and composite (from lamination theory), 
and composite properties represented the actual materials in the areas 
indicated (Fig. 30). The tapered rail configuration was represented by 


^ Personal conversation with Dr. Whitney indicated that his solution 
appeared to be unstable for the [^45]^ laminate. 



FIGURE 30. AXIALLY LOADED RAIL SHEAR SPECIMEN 
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FIGURE 31, TAPERED RAIL SHEAR SPECIMEN 
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512 NODES CD STEEL 



FIGURE 32. OFFSET LOADED RAIL SHEAR SPECIMEN 
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the same mesh as the nontapered rail configuration and the rail thick- 
ness, T^, was tapered to zero at the free edge of the specimen (Fig. 

31). Material properties for the varying rail thickness with uniform 
composite thickness laminate were obtained from lamination theory at 
appropriate intervals along the rails. The finite element representa- 
tion for the shear specimen with offset loading is shown in Fig. 32, 

The dimensions C, F2, T^, and T^ are the same as those of the nontapered 
rail (Fig. 30) and in addition F^ = M. The geometry and loading of this 
specimen are similar to the rail shear specimens used by Whitney et al 
[ 8 ]. 

The normalized shear stress distributions for the three rail configur- 
ations with a [0] laminate (Figs. 33c, 34c, and 35c) indicate a uniform 
shear stress over most of the specimen. The nontapered rail configura- 
tion produces a value of normalized closest to 1.0, at the center of 
the specimen. In all cases the maximum shear stress (Table 5) did not 
occur at the center of the specimen. 

The normal stress contour plots indicate that the highest stress 
concentrations occur in the vicinity of the free edge AB for all rail 
configurations studied. The stresses are maximum at the corners A 
and B. Because of the boundary condition = 0 along the free edge AB, 
the maximum value of is some distance in from the corner A along the 
rai 1 -specimen interface AC for the nontapered and tapered rail configur- 
ations. For the offset loaded rail configuration is very low through- 
out the specimen, and, unlike the axially loaded rail configurations, it 
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is a maximum away from the rail-specimen interface AC. 

A characteristic of all rail shear specimens with elastic rails is 
the existence of nonzero normal stresses in the central portion of 
the specimen. These stresses must be nonzero in this region as a 
result of the nonzero axial strain present in both the fixture and the 
specimen as a result of the elastic behavior of the fixture. The 
effect of rail stiffness on at the center of the specimen is clearly 
illustrated by comparing normalized values for the nontapered and 
tapered rail configurations (Table 5). Since the tapered rail stiffness 
at the center of the specimen is half that of the nontapered rail, 
approximately twice the axial load (and hence stress) is transferred 
into the specimen with the tapered rails. 

7.3.2 Thermal Load Effects 

The thermal stresses in a specimen with a [0] graphite-polyimide 
laminate were determined for a temperature change AT = -600°F using the 
finite element mesh shown in Fig. 30 (with F = 0). Rail configurations 
studied included the nontapered and tapered rails described in section 
7.3.1. The stress contour plots for a^, and normalized with re- 
spect to their corresponding ultimate stresses ^^y^ult’ 

(x^y)uit (Appendix C) are shown in Figs. 36 and 37. 

This linear elastic analysis predicts stress concentrations at the 
corners A and B with both and much greater than their respective 
ultimate values. While the magnitude of shear stress is zero in the 
center of the nontapered rail specimen, it is greater than zero in the 
tapered rail specimen. The nonzero shear stress in the center of the 
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tapered rail specimen is attributed to the coarse finite element grid 
since the shear stress must be zero there to satisfy symmetry condi- 
tions. In general the normalized stresses are considerably lower 
than and however, is significant at the center of the specimen 

7.3.3 Other Mechanically Loaded Laminates 

Normalized stress contour plots for [90] and [±45]^ graphite- 
polyimide laminates using the tapered rail configuration are illustrated 
in Figs. 38 and 39, respectively. The finite element grid described in 
section 7.3.1 was used for these results. 

7. 3. 3.1 The [90] Laminate 

The normalized shear stress distribution for the [90] laminate 
(Fig. 38c) is nearly ideal in that the normalized shear stress reaches 
a value of -1.0 very close to the free edge AB and remains very uniform 
over the entire test section. However, the normal stresses are non- 
zero in the test section and high stress concentrations are present in 
the vicinity of corners A and B. 

7. 3. 3. 2 The [±45]^ Laminate 

The shear stress distribution for the [^45]^ laminate is shown in 
Fig. 39c. As was the case for the rigid rail results (section 7.2), the 
tapered elastic rail solution indicates that the magnitude of the nor- 
malized shear stress at the center of the specimen is indeed close to 
-1.0 (Table 5). However, unlike all previous finite element results for 
both rigid and elastic rails, the maximum normalized shear stress was 
located at corner A where the value was -1.37. The contour plots of the 
normal stresses depict relatively low (but significant) values in the 
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test section and maximum values in the vicinity of the corners. The 
distributions of a and a were similar except that a was generally 

X y jr 

lower in magnitude. 

7.4 Summary 

The rail configurations investigated produced regions of uniform 
shear stress over most of the test section for all laminates studied. 
However, this area of shear stress was often accompanied by significant 
normal stresses depending upon the stiffness of the rails, the method of 
load application, and the laminate properties. Significant axial stress 
was always present in the test section for elastic rails axially loaded. 
The magnitude of this axial stress was approximately inversely propor- 
tional to the stiffness of the rails. The offset loaded rail configura- 
tion minimized undesireable normal stresses in the center of the 
specimen. Stress concentrations in the vicinity of the corners were 
found in all specimens. Although tapered rails reduce stresses caused 
by thermal loads, this linear elastic analysis indicates that a non- 
linear analysis is required for a more accurate solution for a 600° 
temperature change. 
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8. CONCLUSIONS 

A linear elastic finite element analysis has shown that the cross 
beam, losipescu, and rail shear specimens produce a region of shear 
stress in the test section. However, effects attributed to laminate 
properties, core material, and rail configurations significantly in- 
fluenced the magnitude and the purity of the shear stress at the center 
of the specimens. In addition, all specimens exhibited stress concen- 
trations in the vicinity of geometric discontinuities. 

The rail shear specimen produced a region of uniform shear stress 
that was not significantly influenced by laminate properties or rail 
configurations. The axially loaded elastic rails resulted in signifi- 
cant normal stresses at the center of the specimen while the offset 
loaded elastic rail configuration substantially reduced these stresses. 
Significant stress concentrations in the vicinity of the free edges were 
found for all laminates and rail configurations; this was particularly 
true for the thermally loaded specimens. 

The losipescu specimen produced uniform, pure shear at the central 
portion of the specimen for all laminates. Limited studies have shown 
that the stress concentrations at the notcn tips and specimen corners 
can be substantially reduced by rounding the notch tips and using 
doublers in the specimen corners. The effects of elastic fixtures and 
thermal loads were minimal in the test section of this specimen. How- 
ever, the magnitude of shear stress in the test section of the losipescu 
specimen was dependent upon the properties of the laminate. Although 
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this specimen is relatively new (and unproven for composites), experience 
from Its use with metals in addition to the analyses in this study indi- 
cate that the losipescu specimen has potential for future use in testing 
of composite materials. 

The cross beam shear specimen was shown to be a viable specimen for 
the [0/90]^ laminate. However, the [±45]^ laminate produced a non- 
uniform shear stress distribution in the test section and this distri- 
bution was influenced by the stiffness of the core material. A pure, 
uniform area of shear stress was not found in the slotted coupon. 
Therefore, it is not recorrenended for determining shear response of 
composites. 

In summary, the linear elastic finite element results have indi- 
cated that several factors influence the stress distribution present in 
the test section of each specimen. The rail shear and losipescu speci- 
mens, despite specific disadvantages, were found to be the most suitable 
specimens for in-plane shear response of composite laminates. 


BIBLIOGRAPHY 


1 . 


2 . 


3. 


4. 


5. 


6 . 


7. 


8 . 


9. 


10 . 


11 . 


R. W. Fenn, Jr. and R. B. Clapper, "Evaluation of Test Variables 
in the Determination of Shear Strength," Proceedings, ASTM, Vol . 
56, Philadelphia, 1956, p. 842. 

W. W. Breindel, C. L. Seale, and R, L. Carlson, "Evaluation of a 
Single-Shear Specimen for Sheet Material," Proceedings, ASTM. 
Vol. 58, 1958, p. 862. ^ 


N. losipescu, "New Accurate Procedure for Single Shear Testing of 
Metals," J. of Materials . ASTM, Vol. 2, Number 3, Sept. 1967,' 
p. 537. 


R. A. Elkin, 6. Fust, and D. P. Hanley, "Characterization of 
Graphite Fiber/Resin Matrix Composites," Composite Materials: 
Testing and Design . ASTM STP 460 . 19S9, 

C. W. Bert, "Experimental Characterization of Composites," 
Composite Materials (edited by L. J. Broutman and R. H. Krock) , 
Vol. 8, 1974, p. 95. 

P. H. Petit, "A Simplified Method of Determining the Inplane 
Shear Stress-Strain Response of Unidirectional Composites," 
Composite Materials: Testing and Design . ASTM STP 460 . 1969, 


M. E. Waddoups, "Characterization and Design of Composites Mater- 
ials," Composite Materials Workshop (edited by S. W. Tsai, J. C. 
Halpin, and N. J. Pagano), Technomic Publishing Co., Stamford, 
Ct., 1968, p. 254. 

J. M. Whitney, D. L. Stansbarger, and H. B. Howell, "Analysis of 
the Rail Shear Test - Applications and Limitations," J. of 
Composite Materials. Vol. 5, Jan. 1971, p. 24. 

L. Rutherford, C. Bossier, and B. Swain, "Shear Properties of 
Boron/Epoxy Composites," General Precision Aerospace Research 
Center, Systems, Inc., Final Report Contract No. 9-69929 to 
Grumman Aircraft Engineering Corp., Bethpage, New York (1968), 

Advanced Composites Design Guide . Structures Division, Air Force 
Flight Dynamics Laboratory, Wright-Patterson Air Force Base, 

Ohio, section 4.2, Test Methods, December, 1975. 

D. F. Sims, "In-Plane Shear Stress-Strain Response of Unidirec- 
tional Composite Materials," J. of Composite Materials. Vol. 7. 
January 1973, p. 124. 


91 


92 


12. E. L. Bryan, "Photoelastic Evaluation of the Panel Shear Test for 
Plywood," Symposium on Shear and Torsion Testing , ASTM STP 289 , 
1961, p. 90. 

13. S. Dastin, G. Lubin, J. Munyah, and A. Slobodzinski , "Mechanical 

Properties and Test Techniques for Reinforced Plastic Laminates," 
Composite Materials: Testing and Design , ASTM STP 460 , Feb., 

1969, p. 14. 

14. M. W. Hyer and D. 0, Douglas, "A Comparison Between Experiment and 
Theory for a Borsic-Aluminum Picture Frame Shear Test," Technical 
Report 76-T18, School of Engineering, Old Dominion University, 
Norfolk, Virginia, Oct., 1976. 

15. R. N. Hadcock and J. B. Whiteside, "Special Problems Associated 
with Boron-Epoxy Mechanical Test Specimens," Composite Materials : 
Testing and Design , ASTM STP 460 , 1969, p. 27. 

16. D. 0, Douglas, D. E. Holzmacher, Z. C. Lane, and E. A. Thornton, 
"Studies in Finite Element Analysis of Composite Material Struc- 
tures," Technical Report 75-M3, School of Engineering, Old 
Dominion University, Norfolk, Virginia, June, 1975. 

17. D. F. Adams and R. L. Thomas, "Test Methods for the Determination 
of Unidirectional Composite Shear Properties," Advances in Struc- 
tural Composites, 12th National SAMPE Symposium, Vol . 12, paper 
AC-5, 1967. 

18. B. W. Rosen, "A Simple Procedure for Experimental Determination of 
the Longitudinal Shear Modulus of Unidirectional Composites," 

of Composite Materials , Vol. 6, October 1972, p. 552. 

19. H. T. Hahn, "A Note on Determination of the Shear Stress-Strain 
Response of Unidirectional Composites," J. of Composite Materials, 
Vol. 7, July 1973, p. 383. 

20. R. R. Rizzo, "More on the Influence of End Constraints on Off- 
Axis Tensile Tests," J. of Composite Materials , Vol. 3, April 
1969, p. 202. 

21. G, L. Richards, T. P. Airhart, and J. E. Ashton, "Off-Axis Tensile 
Coupon Testing," J. of Composite Materials , Vol. 3, July 1969, 

p. 586. 

22. N. J. Pagano and J. C. Halpin, "Influence of End Constraint in 
the Testing of Anisotropic Bodies," J. of Composite Materials , 

Vol. 2, January 1968, p. 18. 


93 


23. R. B. Pipes and B. W. Cole, "On the Off-Axis Strength Test for 
Anisotropic Materials," J. of Composite Materials. Vol . 7. April 
1973, p. 246, 

24. E. M. Wu and R. L. Thomas, "Off-Axis Test of a Composite," Note 
to the Editor, 0 . of Composite Materials . Vol. 2, October 1968, 
p. 523. 

25. C. C. Chamis and J. H. Sinclair, "10° Off-axis Tensile Test for 
Intralaminar Shear Characterization of Fiber Composites." NASA 
TN D-8215, April, 1976. 

26. C. W. Bert, "Static Testing Techniques for Filament-Wound 
Composite Materials," Composites . January 1974, p. 20. 

27. R. R. Rizzo and A. A. Vicario, "A Finite Element Analysis for 

Stress Distribution in Gripped Tubular Specimens," Composite 
Mater ials: Testing and Design (Second Conference!. ASTM STP 

497, 1972. p. 68. ^ 

28. N. L. Hancox, "The Use of a Torsion Machine to Measure the Shear 
Strength and Modulus of Unidirectional Carbon Fibre Reinforced 
Plastic Composites," J. of Materials Science 7 . 1972, p. 1030. 

29. L. B. Greszczuk, "Shear Modulus Determination of Isotropic and 

Composite Materials," Composite Materials; Testino and Desion. 
ASTM STP 460 . 1969, p. 140. — 

30. S. A. Sattar and D. H, Kellogg, "The Effect of Geometry on the 

Mode of Failure of Composites in the Short Beam Shear Test," 
Composite Materials: Testing and Design . ASTM STP 460 . 1969, p. 


31. C. A. Berg, J. Tirosh, and M. Israeli, "Analysis of Short Beam 

Bending of Fiber Reinforced Composites," Composite Materials: 
Testing and Design . ASTM STP 497 . 1972, p7 206. ~~ 

32. R. M. Jones, Mech anics of Composite Materials, Scripta Book Co.. 

Washington, D, C., 1975. ~ “ 

33. C. S. Desai and J. F. Abel, Introduction to the Finite Element 
l^ethod . Van Nostrand Reinholdl Co. , New York, 1972. 

34. 0. C. Zienkiewicz, The Finite Element Method in Engineering 
Science . McGraw-Hill , London , 1971, 

35. E, L. Wilson, K. Bath, and W. P. Doherty, "Direct Solution of 
Large Systems of Linear Equations." Computers and Structures. 

Vol. 4, p. 363, 1974. 


APPENDIX A. 
COMPUTER PROGRAMS 


94 


A. COMPUTER PROGRAMS 


A.l Introduction 

The present study required three distinct FORTRAN programs to 
obtain the levels of efficiency and analysis desired. The finite element 
program described in section A. 2 is the basic program. This program 
performs the finite element solution independently of the other pro- 
grams. However, for an accurate finite element solution, many intrin- 
sically repetitive node and element cards are required which collectively 
represent the geometry of the structure. Since repetitive numerical 
procedures are simplified when programmed, a mesh generation program was 
developed for the requirements, of the study as described in section A. 3. 

The number of elements in a finite element representation produce 
a proportional number of stress values which are provided in the form of 
tabular output by the basic finite element programs. The program 
described in section A. 4 conveniently reduces this data by contour 
plotting the tabulated stress values. 

In practice the finite element mesh is generated first and the 
resulting node and element card images are written on an indirect access 
file. Since the node and element card images constitute most of the 
input to the finite element program (see section A. 2) the file is 
interactively edited to become an input data file for the finite element 
program. Similarly, the tabulated stress values from the finite element 
program may be stored on an indirect access file. This file constitutes 
a large portion of the input for the stress contour plotting program 
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(see section A. 4) and it also is edited to become an input file to the 
contour plotting program. 

A. 2 The Finite Element Program 

A. 2.1 General 

The finite element method is a powerful tool for approximating the 
state of stress in virtually any structure subjected to loads. A re- 
quirement for utilizing such a method is the availability of a high 
speed digital computer and a program to carry out the solution procedure. 
This section describes such a FORTRAN program. 

The use of the finite element method for solutions to structural 
problems has been well established by Desai and Abel [33], Zienkiewicz 
[34], and others. This particular program is an extensively modified 
version of the FORTRAN code by Desai and Abel. The program uses the 
displacement method for solution (see section 3.2). Both constant 
strain triangular elements and quadrilateral elements composed of four 
constant strain triangles are available for either plane stress or plane 
strain analysis. Capability added to the program by the author includes 
orthotropic material properties, orthotropic thermal effects, and the 
implementation of an equation solving subroutine developed by Wilson et 
al [35]. The improved equation solver allows the program to solve a 
very large system of linear equations in blocks residing on slow speed 
storage, as opposed to the more restricted, completely in core solution 
technique. 

The job flow is one pass since linear elastic problems are directly 
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solved. After material property, node, element, and load information is 
input, a stiffness matrix for each element is computed. All elemental 
stiffness matrices are then assembled into a global stiffness matrix, 
and the load vector is initialized. The unknown displacement vector is 
determined using Gauss elimination techniques on only non-zero stiffness 
terms. Once the global displacements are computed, elemental strains 
are determined. The elemental stresses are then computed from the 
elemental constitutive matrices and the elemental strains. 

A. 2.2 Input Description 
Column Contents 

Title card (I5,9A8) 

1-5 NPROB = Problem number. If NPROB > 100, 

stresses are written on unit 2, 
This option is used with the 
stress contour plotter (see 
section A. 4). If NPROB = 0, 
no more problems are run. 


6-77 

A(I) 

= Title information 

Control card (415) 

1-5 

NNP 

= Number of nodes 

6-10 

NEL 

= Number of elements 

11-15 

NMAT 

= Number of materials 

16-20 

NOPT 

= 1 for plane strain analysis 
= 2 for plane stress analysis 
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Laminate material property cards {7E10.0) 
(for material I) 


1-10 

11-20 

21-30 

31-40 

41-50 

51-60 

61-70 


6X(I) 

= Modulus of elasticity in X 
direction 

EY(I) 

= Modulus of elasticity in Y 
direction. IF (EY(I) .EQ.O. ) 
EY(I) = EX(I) 

PRXY(I) 

= Poisson's ratio, v . IF 

xy 

(PRXY(l).EQ.0.)PRXY(I) = EX(I)/ 
2./GXY(r)-l. 

GXY(I) 

= Modulus of rigidity for materi- 
al I. IF (GXY(I).EO.O.) 

GXY(I) = EX(l)/(l.+PRXY(I))/2. 

TH(r) 

= Thickness for material I. 


ALPHA(1,I)= Coefficient of thermal expansion 
in X direction for material I. 
ALPHA(2,I)= Coefficient of thermal expansion 
in Y direction for material I. 


Nodal location and mechanical load cards (2I5,4E10.0) 

1-5 N = node number 

6-10 KODE(N) = 0 for ULX(N) and VLY(N) pre- 

scribed loads 

= 1 for ULX(N) prescribed displace- 
ment and VLY(N) prescribed load. 
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= 2 for ULX(N) prescribed load 
and VLY(N) prescribed displace- 
ment 


= 3 for ULX(N) and VLY(N) pre- 
scribed displacements 


11-20 

X(N) = X coordinate 

21-30 

Y(N) = Y coordinate 

31-40 

ULX(N) = X displacement or load 

41-50 

VLY(N) = Y displacement or load 

Element and thermal 

load cards(6I5.E10.0) 

1-5 

M = Element number 

6-10 

IE(M,1) = First node of element 

11-15 

IE(M,2) = Second node of element 

16-20 

IE(M,3) = Third node of element 

21-25 

IE(M,4) = Fourth node of element. For 

triangular elements, IE(M,3) = 
IE(M,4) 

26-30 

IE(M,5) = Material properties associated 

with element M 

31-40 

TEMP(M) = Change in temperature of 


element after displacement 
boundary conditions are imposed. 
IF (TEMP{M).EQ,0.)TEMP(M) = 
TEMP(M-l) 
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Comments 

Node and element cards must be in ascending consecutive order. 

Nodes associated with triangular elements that are not ordered counter- 
clockwise cause a warning message, but are reordered counterclockwise 
internally. Quadrilateral elements must be ordered counterclockwise. 

If the mesh generation program (section A. 3) is used to supply node and 
element cards, these cards (or indirect access file) must be altered to 
include boundary condition (KODE), mechanical load (ULX(N),VLY(N)), 
elemental material (IE(M,5)), and elemental temperature (TEMP(M)) 
information. To terminate execution, a blank card is inserted after the 
last element and thermal load card of the last model. If NPROB is not 
zero the program expects to read a complete data set for another model. 

A. 3 The Mesh Generation Program 

A. 3.1 General 

This program generates groups of elements and nodes. In a group, 
all elements are either identical to one another or have repetitive 
geometric characteristics. Many groups of elements make up the total 
model. A group of elements may be located anywhere in X-Y coordinates 
by several methods. The dimensions, orientation, number of elements, 
and type of elements for each group are specified by the user. Data for 
each element group is read sequentially so that data to generate speci- 
fic element groups may be added or deleted from the entire data set. 
Specified groups of elements may be scaled, rotated, or repositioned so 
that element groups frequently needed in different meshes need only be 


101 


checked once and can be inserted as needed into different meshes. 

Input data to the mesh generator is comprised of four categories. 
The first category of input is the generation initialization and termin- 
ation data. The initialization data is at the beginning of the data set 
and is used only once for the purpose of initializing title information 
and reference frames. At the end of the data is the card indicating 
that the mesh generation commands are complete so that the program can 
perform its final computations. 

All other categories of input commands may be used anywhere at any 
time between the initialization and termination data. The second 
category of data consists of element group generation data sets. These 
data sets generate specific groups of elements. An element group gener- 
ation data set can control every parameter associated with the genera- 
tion of that particular group, or certain parameters may be passed to it 
from the proceeding group. The third category of input alters the 
reference frame of one, several, or all groups that constitute the 
complete model. General capabilities include translation, rotation, and 
scaling of any desired element groups. The last category of input data 
has the capacity to reposition nodes. While other input commands 
operate by element groups, the nodal repositioning commands operate on 
all nodes that have been previously generated. The only capacity for 
nodal repositioning in the present version is that of repositioning 
nodes to lie in a precise circular arc. This avoids the necessity of 
manually positioning individual nodes that will ultimately lie on some 
circular arc in the complete model. 
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A. 3. 2 Generation Initialization and Termination Data Input 
Column Contents 

Title card (10A8) 

1-80 CARD(l) = Title information 

Coordinate Initialization Card (6F10.0»2I5) 


1-10 

ORIGIN(l) 

= Absolute X coordinate of initial 
global and local coordinate 
system origin 

11-20 

0RIGIN(2) 

= Absolute Y coordinate or initial 
global and local coordinate 
system origin 

21-30 

SCALE(l) 

= X scaling dimension for internal 
node numbering 

31-40 

SCALE(2) 

= Y scaling dimension for internal 
node numbering 

41-50 

FACT(l) 

= X global coordinate scale factor 

51-60 

FACT(2) 

= Y global coordinate scale factor 

61-65 

MSCAN 

= Node duplication elimination 
number 

66-70 

NSV 

= 0 for numbering nodes by as- 


cending X then Y location 
= 1 for numbering nodes by as- 
cending Y the X location 
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Comments 

The ORIGIN(I) variables initialize the origin of the global coordin- 
ates in terms of absolute (final) coordinates. Since the ORIGIN(I) 
fields have the same purpose as the global origin reset cards, the ex- 
act use of these fields are described in section A. 3. 4. 

The node numbers are internally numbered by absolute X and Y loca- 
tion. Such node numbers are determined by the equations 

IX = X/SCALE(l)*10000.-.5 
lY = Y/SCALE(2)*10000.-.5 
IF(NSV.EQ.0)N0DE = IY*10000+IX 
IF{NSV.EQ.1)N0DE = IX*10000+IY 

where X,Y are node locations. 

NODE is the internal node number. The nodes are internally numbered 
in this manner for two reasons. Due to roundoff errors caused by 
reference frame alterations, a previously generated node does not neces- 
sarily have the exact coordinate location of another node intended to 
have the same location. For this reason, duplicate node positions are 
not eliminated by a check for X-Y absolute coordinate duplication. 
SCALE(I) must be slightly larger than the expected maximum X-Y node 
coordinates. NODE is then determined for each node, and any node whose 
internal IX or lY is equal to or one different from another IX or lY 
is eliminated. After all nodes and elements have been generated, 
the NODE numbers are put in ascending order and renumbered 1 ,2,3, . . . ,NNP, 
where NNP is the total number of node points. This procedure minimizes 


104 


the bandwidth of many meshes since nodes are renumbered by ascending 
X coordinate then ascending Y coordinate locations (for NSV=0). 

The FACT(I) variables initialize the coordinate multiplication 
factors for all dimensions with respect tc any global coordinate system. 
Since the FACT(I) fields have the same purpose as the coordinate scale 
reset cards, the exact use of these fields are described in section A. 3. 4. 

The integer MSCAN is used to minimize scanning procedures for 
duplicate NODE numbers as element groups are processed. For example, 
if an element group with ten boundary nodes is being added to an exist- 
ing model whose last ten generated nodes will coincide with these 
boundary nodes, MSCAN need only be ten. However, if one hundred new 
nodes are generated between the boundary nodes that have the same loca- 
tion, MSCAN should be at least one hundred and ten. Since all duplicate 
nodes are eliminated at the time the model is completed, MSCAN can be 
zero for any model. However, duplicate node information tends to fill 
up storage arrays rather quickly so the value of this variable should 
always be overestimated to minimize storage requirements. 

Column Contents 

Termination card (15) 

1-5 KIND =5. This stops the program from 

reading additional group genera- 
tion, reference frame altera- 
tion, and nodal repositioning 
cards. At this point the pro- 
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gram initiates final processing 
of the completed mesh. This 
card may be last in the input 
deck (or file) or placed be- 
tween any of the repeatable 
data sets for mesh debugging 
purposes 

A. 3. 3 Element Generation Input 
Type 1 Rectangular Element Generation Cards 
Card 1 (I5,9A8) 

KIND = 1 for group type 1, rectangular 
element generation input to 
f ol 1 ow 

CARD(I) = Title for element group 

Card 2 (3I5,F10.0) 

^■5 M = The number of rectangular ele- 

ments in the X direction, see 
figure A-la 

6-10 N = The number of rectangular ele- 

ments in the Y direction 

11-15 IRC = The reset code for A(l) and A(2) 

16-25 THETA = The angle of rotation about the 

local origin 

Card 3 (2F10.0,I5) 
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1-10 

11-20 

21-25 


A(1) or = Total reset or multiplying 

AX factor of A(1 ) 

A(2) or = Total reset or multiplying 

AY factor of A(2) 

JRC = Reset code for XY(1) and 

XY(2), the local origin loca- 
tion in terms of global coor- 
dinates 


Card 4 (2F10.0) 


1-10 

XY(1) or 

= Total or incremental X local 


TEMPO ) 

origin location in terms of 



global coordinates 

11-20 

XY(2) or 

= Total or incremental Y local 


TEMP{2) 

origin location in terms of 



global coordinates 

Comments 



For the A(I) 

dimension reset option 

IRC, the following numbers are 

admissable 

IRC 

= 0 for no reset, A(l) and A{2) 


are used from the previous 
element group 


= 1 for a scalar reset of A(l), 
AX is read, A(1) = AX*A(1)* 
FACTO ) 
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= 2 for a scalar reset of A(2), 
AY Is read, A(2) = AY*A(2)* 
FACT{2) 

= 3 to reset both cases above 
= -1 for a total reset of A(l), 
A(1 ) is read, A(1 ) = A(1 )* 
FACT(l) 

~ -2 for a tatal reset of A(2),,_ 
A(2), is read, A(2) = 
A(2)*FACT(2) 

= -3 for a total reset of A(l) 
and A(2) 

The program omits fields that are not applicable following resets. If 
IRC = 0, the following field on card 3 becomes (15). If IRC 1, the 
following field on card 3 becomes (FIO.0,15), and similarly for other 
cases. 

For the local origin reset the following options are admissable 

JRC = 0 for the default change in 

local origin with respect to 
the global coordinate system 
= 1 for an incremental reset of 
XY(1), TEMP(l) is read, XY(1) = 
XY(1)+TEMP(1)*FACT(1)+0RIGIN(1) 
= 2 for an incremental reset of 
XY(2), TEMP(2) is read. XY(2) = 
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la ' g 


M=3 , N= 2 


■FINAL LOCAL ORIGIN 


ORIGIN (I) 


XY{2) 


_1 

A(2) 

"T 


ORIGIN (2) 




INITIAL LOCAL ORIGIN - 


a. GROUP I - RECTANGLES 


a-ABSOLUTE 
g? GLOBAL 



KIND = 2 


b. GROUP 2 - ISOSCELES TRIANGLES 


XY(I)^ 


XY (I)i 



KIND = - 3 


KIND = 3 


c. GROUP 3- TRANSITION "RECTANGLES" 


TEMP(I) 1 


XY(I) 


KIND : - 4 


KIND= 4 


d. GROUP 4- IRREGULAR ELEMENTS 


FIGURE A-1 ELEMENT GROUPS 


XY(2)+TEMP(2)*FACT(2)+0RIGIN(2) 
= 3 to reset both cases above 
= -1 for a total reset of XY(1)> 
XY(1) is read, XY(1) = XY(1)* 
FACT(1)+0RIGIN{1) 

= -2 for a total reset of XY(2), 
XY{2) is read, XY(2) = XY(2)* 
FACT(2)+0RIGIN{2} where XY(I) 
become absolute coordinates. 

= -3 for total reset of XY(1) 
and XY(2) 

All element group generation cards may be used anywhere and any number 
of times after the initialization cards. 


Type 2 Isosceles Triangle Element Generation Cards 
Card (I5,9A8) 

kind = 2 or -2, depending on the 
generation mode desired for 
isosceles triangle input, see 
figure A-lb 

CARD(I) = Title information for element 
group 


Card 2 (3I5,F10.0) 
1-5 


= The number of isosceles triangles 


in the X direction 


no 


6-10 

n-15 

16-25 

Card 3 (2F10.0,I5) 
1-10 

11-20 

21-25 

Card 4 (2F10.0) 
1-10 

11-20 


Comments 


N = The number of isosceles tri- 

angles in the Y direction 
IRC = The reset code for A(I), see 

type 1 description 

THETA = The angle of rotation about 
the local origin 


A(l) or 

= Total reset or multiplying 

AX 

factor of A(1 ) 

A(2) or 

= Total reset or multiplying 

AY 

factor of A(2) 

JRC 

= Reset code for XY(I), the 


absolute coordinate location 


for the local origin, see 


group 1 


XY(1) or 

= Total 

or incremental 

X local 

TEMP(1 ) 

origin 

reset value 


XY(2) or 

= Total 

or incremental 

Y local 

TEMP(2) 

origin 

reset value, 

see 


group 

1 



See figure A-lb for isosceles triangle specifications. The 




in 



description of cards 3 and 4 

are identical to those of element group 


1. 




Type 3 Transition "Rectangle' 

Generation 

Cards 


Card 1 (I5,9A8) 




1-5 

KIND 

3 or -3, depending on the mode 




for transition rectangles, see 

i 



figure A-lc 

1 

6-77 

CARD(I) 

Title information for element 

1 

i 

1 



group 

1 

Card 2 (2I5,F10.0) 



1 

! 

1-5 

M 

The number of transition 

1 

i 



rectangles in the X direction 

! 

6-10 

IRC 

The reset code for A(I) 


11-20 

THETA 

The angle of rotation about 

1 



the local origin 

I 

Card 3 (2F10.0,I5) 




1-10 

A(1) or = 

Total reset or multiplying 

! 


AX 

factor of A(1 ) 

; 9 

11-20 

A(2) or 

Total reset or multiplying 

«4 


AY 

factor of A(2) 


21-25 

JRC 

Reset code for XY(1) and XY(2), 




the local origin location 


Card 4 

) 

1 

I 

/ 
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1-10 

XY{1) or = 

Total or incremental X local 


TEMP(l) 

origin reset value 

11-20 

XY(2) or “ 

Total or incremental Y local 


TEMP(2) 

origin reset value 

Type 4 Irregular 

Element Generation Cards 


Card 1 (I5,9A8) 



1-5 

KIND 

4 for irregular triangle gen- 



eration 


= 

-4 for irregular quadrilateral 



generation 

6-77 

CARD(I) = 

Title information for element 



group 

Card 2 (215) 



1-5 

NELODD 

The number of irregular ele- 



ments to be generated, see 



figure A-ld 

6-10 

ORC 

The local origin reset code 

Card 3 (2F10.0) 



1-10 

XY(1) or = 

Total or incremental X local 


TEMP(l) 

origin reset value 

11-20 

XY(2) or = 

Total or incremental Y local 


TEMP{2) 

origin reset value 
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# 


Card 4 (6F10.0) 


1-10 

TEMP(l) 

11-20 

TEMP{2) 

21-30 

TEMP(3) 

31-40 

TEMP(4) 

41-50 

TEMP(5) 

51-60 

TEMP(6) 


= Change in X coordinate from 
local origin for second node 
on element 

= Change in Y coordinate from 
local origin for second node 
on element 

= Same as TEMP(l) for third 
node 

= Same as TEMP(2) for third 
node 

= Same as TEMP(l) for fourth 
node 

= Same as TEMP(2) for fourth 
node 


Card 5 (2F10.0) 

1-10 TEMP(l) 


11-20 


TEMP(2) 


A. 3.4 Reference Frame Alteration 


= Change in X coordinate of the 
local origin after each irregu 
lar element generation 
= Change in Y coordinate of the 
local origin after each irregu 
lar element generation 
Data Sets 


Global Origin Reset Cards 
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Card 1 (I5,9A8) 

1-5 
6-77 

Card 2 (2F10.0) 

1-10 


11-20 


Comments 

These cards may be inserted anywhere after the initialization cards 
and may be used any number of times. The local coordinate position is 
always initially set to the ORIGIN(I) position. All local origin 
changes and resets are done with respect to the global coordinate system. 

Global Coordinate Scale Factor Reset Cards 

Column Contents 

Card 1 (I5,9A8) 

1-5 KIND = 6 for global scale factor reset 


KIND = -1 for global origin reset 

CARD(I) = Title information 

ORIGIN(l) = Location of X global origin in 
absolute coordinates and ini- 
tial location of the local 
coordinate system 

0RIGIN(2) = Location of Y global origin in 
absolute coordinates and ini- 
tial location of the local 
coordinate system 


6-77 


CARD(I) = Title information for scaled 
element groups 
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FACT(l) = Scale factor for global coor- 
dinate in X direction 
FACT{2) = Scale factor for global coor- 
dinate in Y direction 

Comments 

FACT(I) scales all A(I), AX, AY, XY(I), and TEMP(I), input from 
the element group generation data sets (see section A. 3.3). 

Input for Rotating Several Element Groups 

Column Contents 

Card 1 (I5.9A8) 

1-5 KIND = 8 for element groups following 

to be rotated 

6-77 CARD(I) = Title information 

Card 2 (5F10.0) 

1-10 XTH = Rotation in degrees about the 

X axis at XPIVOT and YPIVOT 

11-20 YTH = Rotation in degrees about the 

Y axis at XPIVOT and YPIVOT 

21-30 ZTH = Rotation in degrees about the 

Z axis at XPIVOT and YPIVOT 

31-40 XPIVOT = Global X location for the 


Card 2 (2F10.0) 
1-10 

11-20 
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point about which nodes are 
rotated 

41-50 YPIVOT = Global Y location for the 

point about which all nodes 
are rotated 

Rotated (15) 

= -8, to rotate nodes of element 
groups generated between card 
2 and this card 

Comments 

XPIVOT and YPIVOT are not scaled by FACT(I). Rotation commands may 
not be nested and for every KIND = 8 there must be a corresponding KIND 
= -8. All rotations executed with these commands are carried out 
independently of rotations of individual element groups. 

A. 3. 5 Node Repositioning Data Sets 
Column Contents 

Card 1 (I5,9A8) 

1-5 KIND = 7 for node repositioning data 

to follow 

6-77 CARD(l) = Title information 

Card 2 (7F10.0) 


Card to Terminate Element Groups to be 
1-5 KIND 


1-10 


XO 


= Absolute X location for the 
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11-20 


21-30 


31-40 

41-50 


51-60 


61-70 


center of the circular segment 
on which' repositioned nodes 
are to lie 

YO = Absolute Y location for the 

center of the circular segment 
on which repositioned nodes 
are to lie 

THI = Initial angle from the X axis 

for circular segment 
THF = Final angle from the X axis 

for circular segment, THF > THI 
TDEG = Tolerance to include nodes lo- 
cated near angles THI and THF. 
Nodes that are repositioned lie 
between the angles THI-TDEG and 
THF+TDEG 

RAD = Radius of the circular segment 

about XO, YO to which nodes 
will be repositioned 
TOL = Tolerance to include nodes 

located between a radius of 
RAD-TOL and RAD+TOL to be 
repositioned along the radius 


RAD 
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A. 4 The Plotting Program 
A. 4.1 General 

In order to efficiently evaluate finite element meshes and stress 
fields, it is desirable to have a versatile computer program capable of 
plotting finite element meshes and stress contour plots either on paper 
or at an interactive graphics terminal. This section describes such a 
finite element mesh and stress contour plotting computer program. 

The finite element mesh plotting algorithm simply plots the speci- 
fied elements of the mesh to a specified plot width. Nodes, elements, 
or both may be numbered. Numbering schemes are chosen by the user. 

The stress contour plotting algorithm plots the outline of the 

mesh and plots lines of constant o , o , and/or t throughout the 

X y xy 

interior of the mesh outline. The number of stress contours are speci- 
fied by the user. Stress contour values may be input, internally 
generated, normalized by input values of stress, or internally norma- 
lized. Any stress contour values generated may be punched for future 
use. 

A.4.2 Input Description 
Column Contents 

Control Card (315) 


1-5 

NNP 

= Number of nodes 

6-10 

NEL 

= Number of elements 

11-15 

NCASE 

= Number of plot cases 
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Node Cards (10X,2F10.0) 


11-20 

XY(I.l) 

= X coordinate of node I. 

21-30 

XY{I,2) 

= Y coordinate of node I. 

Element Cards (5X,4l5) 



6-10 

IELE(I,1) 

= First node number of element I 

11-15 

IELE(I,2) 

= Second node number of element 
I 

= Third node number of element I 

16-20 

IELE(I,3) 

21-25 

IELE(I,4) 

= Fourth node number of element 
I 

Stress Cards (3E20.0) 



1-20 

SIG(I,1) 

= at the centroid of element 

I 

= Oy at the centroid of element 
I 

= at the centroid of element 

21-30 

SIG(I,2) 

31-40 

SIG(I,3) 


I 


Comments 

There are NNP node cards, node numbers are in ascending order from 
1 to NNP. There are NEL element cards in ascending order from 1 to NEL. 
IELE(I,3) = IELE(I,4) for triangular elements. If only a grid is to be 
plotted, NCASE is negative and the number of plot cases run are 
lABS(NCASE). No stresses are read for a negative NCASE. There are NEL 
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stress cards, one associated with every element. Node and element cards 
(or records) are usually generated with the mesh generator and the 
stress cards are usually generated by the finite element program 
(NPROB > 100, see section A. 2). 


Plot Case Cards 


Column 


Contents 

Title Card (10A8) 

1-80 

TITLE(I) 

= Title for plot case to follow 

Control Card (5F10, 

.0,315,611 ) 


1-10 

XMIN 

= Minimum X dimension for nodes 

11-20 

XMAX 

plotted 

= Maximum X dimension for nodes 

21-30 

YMIN 

plotted 

= Minimum Y dimension for nodes 

31-40 

YMAX 

plotted 

= Maximum Y dimension for nodes 

41-50 

WIDTH 

plotted 

= Width of the plotted mesh 

51-55 

NEON 

= The labeling index for node 


' 

numbering on mesh plots 

56-60 

NEOE 

= The labeling index for element 

61-65 

NCONT 

numbering on mesh plots 
= The number of stress contours 
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to be plotted 

66 IGRID = 0 for no grid plot 

= 1 for grid plot 

67 ISIGX = 0 for no stress contour 

plot 

= 1 for o stress contour plot 
= 2 for stress contour plot 
with all plotted stresses 
normalized by normalization 
elements 

= 3 for stress contour plot 
with all plotted stresses 
normalized by an input value 
of stress 

68 ISIGY = Same options as field 67 (ISIGX) 

for 

69 ITAUXY = Same options as field 67 

(ISIGX) for 

70 ICOVAL = 0 for elementwise stress con- 

tour value generation 
= 1 for input stress contour 
values 

= 2 for stress contour values of 


constant stress increment 
generation 
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71 ICOPUN = 0 for no punched stress contour 

values 

= 1 for punched stress contour 
values 

Normalizing Stress Card (FIO.O) 

1-10 AVE(4) = Stress which all others will 

be normalized by. This card 
is needed only if ISIGX, 

ISIGY or ITAUXY is 3. 

Stress Normalization Element Cards (15, FIO.O, 1215) 

1-5 ICON = 0 for no more elements to be 

added 

= +N for every Nth element to be 
averaged 

= -1 for each following element 
to be averaged 


6-15 

FACT 

= Element's weighting factor 

16-80 

NOPAIR 

= Element numbers 


(1,1) 



Plotted Element Cards (1415) 

1-5 ICON = 0 for no more elements to be 

added 

“ +N for every Nth element to be 


plotted 

= -1 for each following element 
to be plotted 


6-75 

NOPAIR = Element numbers 

(M) 

Stress Value Cards 

(3E20.0) 

1-20 

AVE(I,1) = cj^ stress value of contour I 

21-40 

AVE(I+1,1)= o stress value of contour I+l 

41-60 

AVE(I+2,1)= stress value of contour 1+2, 

repeat cards until NCONT values 
are read 

1-20 

AVE(I,2) = Oy stress value of contour. I 

21-40 

AVE(I+1,2)= Oy stress value of contour I+l 

41-60 

AVE(I+2,2)= Oy stress value of contour 1+2, 
repeat cards until NCONT values 
are read 

1-20 

AVE(I,3) = T stress value of contour I 
xy 

21-40 

AVE(I+1,3)= T stress value of contour I+l 
xy 

41-60 

AVE(I+2,3)= T stress value of contour 1+2 
xy 

repeat cards until NCONT values 
are read. If ISIGX, ISIGY, 
and/or ITAUXY are 0, their 
corresponding stress contour 
values are not needed. Contour 
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stress value cards are only 
read for IC0VAL=1 • 


Comments 

All plot case cards are repeated NCASE times. For automatic 
scaling, XMIN, XMAX, YMIN, and YMAX are zero. The smallest of XMAX-XMIN 
and YMAX-YMIN is plotted across the plotting drum. For nonzero values 
of XMIN, XMAX, YMIN, and YMAX the plot will be distorted accordingly. 
WIDTH is the maximum dimension across the plotting drum for the plotted 
elements when XMIN, XMAX, YMIN, and YMAX are zero. A poor choice of 
XMIN, XMAX, YMIN, and YMAX may cause WIDTH to be significantly different 
from the expected width. 

The value of NEON determines which nodes are numbered according to 
the formula 

number of labeled node = 1+K/NEON (NEONj^O) 
where K = 1,2,3,..., NNP. 

No nodes are numbered if NEON = 0. Element numbering is determined in 
the same manner with NEOE. Element numbers are underlined. 

Elementwise stress contour value calculations provide stress con- 
tours in the critical, refined mesh areas. For ICOVAL = 0 or 2, the 
elemental stresses are first put in ascending order in array X. The 
stress contour values 1 through NCONT-1 are calculated with the equations 


K = NEL/NCONT 
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AVE(I) = (X(I+L)+X(I+l+L)+...+X(I+K+L))/K 
where L = K*(M-1) for M = 1 ,2,3, . . . ,NC0NT-1 . 

The last stress contour value is the average of the stresses at the re- 
maining elements. 

Stress contour values of constant stress increments for ICOVAL = 2 
are generated with the equations 

DX = (X(NEL)-X(l))/(NC0NT+2} 

AVE(I) = X(1)+DX*I for I = 1 ,2,3, . . . ,NC0NT. 


The stress normalization cards determine a normalizing a , o , or t 

a x’ y’ xy 

stress value from the stresses on one or several elements weighted by 
FACT. For example, the records 


1 2. 1 3 6 8 10 

2 1 . 11 21 

-1 .5 17 18 35 42 56 

0 

would cause the (for ITAUXY = 2) stress value of the indicated 

element numbers to be weighted as follows 

A = (o3+o^+a5+ag+a^+ag+ag+a.,Q)/2. 

B = (o] ■i+f’i3+Oi5+OT7+'7-| 9+021 )/l • 

C = (oi7+Oig+035+042‘'"5g)/0.5 

= ^VUm= (A.B. 0 / 18 ., 
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where 


a. are t stresses of the element. 
1 xy 


If neither ISIGX, ISIGY, nor ITAUXY are equal to 2, these cards 


are not included. Element numbers on individual cards must be in 
ascending order. 


APPENDIX B 
PROGRAM LISTING 

B.l THE FINITE ELEMENT PROGRAM 
B.2 THE MESH GENERATION PROGRAM 
B.3 THE PLOTTING PROGRAM 

Listings of these programs are available upon request. 
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APPENDIX C. 
MATERIAL PROPERTIES 
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Product Liability: Industry Study. Volumes I and II 

PB-265 542/PSK 886 p PC$24.00/MF$3.00 

An Atlas of Radiation Histopathology 

TID-26676/PSK 234 p PC$7.60/MF$3.00 

The Medical School Admissions Process. A Review 
of the Literature, 1955-76 

PB-263 962/PSK 193 p PC$9.00/MF$3.00 

Integrated Energy Vocabulary 

PB-259 000/PSK 459 p PC$22.50/MF$22.50 

Nuclear Air Cleaning Handbook. Design, Construc- 
tion, and Testing of High-Efficiency Air Cleaning 
Systems for Nuclear Application 

ERDA-76/21/PSK 300 p PC$9.75/MF$3.00 

Evaluation of the Air-To-Air Heat Pump for 
Residential Space Conditioning 

PB-255 652/PSK 293 p PC$11.00/MF$3.00 

Federal Information Processing Standards Register: 
Guidelines for Documentation of Computer 
Programs and Automated Data Systems. Category: 
Software. Subcategory: Documentation 

FIPS-PUB38/PSK 55 p PC$5,25/MF$3.00 

Life Cycle Costing Emphasizing Energy Conserva- 
tion: Guidelines for Investment Analysis 

ERDA-76/130/PSK 117 p PC$6.50/MF$3.00 

Analysis of Projected Vs. Actual Costs for Nuclear 
and Coal-Fired Power Plants 

FE-2453-2/PSK 38 p PC$4.50/MF$3.00 

Quality Criteria for Water 

PB-263 943/PSK 537 p PC$15.50/MF$3.00 

Passive Solar Heating and Cooling Conference and 
Workshop Proceedings May 18-19, 1976, 

University of New Mexico, Albuquerque, New 
Mexico 

LA-6637-C/PSK 352 p PC$12.50/MF$3,00 

Predicting the Performance of Solar Energy Systems 

AD-A035 608/PSK 45 p PC$4.50/MF$3.00 

NIOSH Analysis Methods for Set J 

PB-263 959/PSK 128 p PC$7.25/MF$3.00 

Flat-Plate Solar Collector Handbook: A Survey of 
Principles, Technical Data and Evaluation Results 

UCID-17086/PSK 96 p PC$6.00/MF$3.00 

Coal Preparation Environmental Engineering Manual 

PB-262 716/PSK 729 p PC$21.50/MF$3.00 


Administration 

Development of Executive Success 
Theodore E. Elsasser. 

Naval Postgraduate School Monterey Calif Sep 76, I87p 
AD-A035 839/OPSK PC$9.00/MF$3.00 

Behavior and Society 

Interagency Task Force on Product Liability. Insurance Study 
McKinscy and Co., Inc., New York. Jan 77, 263p 
PB-263 600/9PSK PCS10.75/MF$3.00 

Biomedical Technology 
and Engineering 

Economical Multifactor Designs for Human Factors 
Engineering Experiments 
Charles W. Simon. 

Hughes Aircraft Co Culver City Calif Engineering Equipment 
Div Jun 73, 190p 

AD-A03S 108/0PSK PC$9.00/MF$3.00 


Building Technology 

Criteria for Selection and Design of Residential Slabs-on- 
Ground 

National Research Council, Washington, D.C. Building 
Research Advisory Board. 1968, 303p 
PB-261 551/6PSK PCS! 1.75/MF$3.00 

Business and Economics 


Interagency Task Force on Product Liability Legal Study 

Research Group, Inc., Charlottesville, Va. Jan 77, 1274p in 
7v 

PB-263 601/7PSK PC$34,00/MF$3.00 


Chemistry 

Research Program on Hazard Priority Ranking of 
Manufactured Chemicals (Chemicals 21-40) 

S. L. Brown, F. Y. Chan, J. F. Jones, D. H. Liu, and K. E. 
McCaleb. 

Stanford Research Inst., Menlo Park, Calif. Apr 75, 195p 
PB-263 162/OPSK PC$9.00/MF$3.00 

Research Program on Hazard Priority Ranking of 
Manufactured Chemicals (Chemicals 41-60) 

S. L. Brown, F. Y. Chan, J. L. Jones, D. H. Liu, and K. E. 
McCaleb. 

Stanford Research Inst., Menlo Park, Calif. Apr 75, 191p 
PB-263 163/8PSK PC$9.00/MF$3.00 

Research Program on Hazard Priority Ranking of 
Manufactured Chemicals (Chemicals 61-79) 

S. L. Brown, F. Y. Chan, J. L. Jones, D. H. Liu, and K. E. 
McCaleb. 

Stanford Research Inst., Menlo Park, Calif. Apr 75, 198p 
PB-263 164/6PSK PC$9.00/MF$3.00 

Civii Engineering 

Asphalt Concrete Overlays of Flexible Pavements. Volume 2. 
Design Procedures 

Austin Research Engineers, Inc., Tex. Jun 75, 96p 
PB-263 433/SPSK PC$6.00/MF$3.00 


Communication 

Problems of Space and Terrestrial Microwave Propagation 
W. Riedler, and W. Lothaller. 

European Space Agency, Paris (France). May 76, 286p 
Partly in English and French. 

N77-132.38/9PSK PCSI 1 .00/MF$3.00 

Modulation Studies for Direct Satellite Communication of 
Voice Signals 
Hiroshi Akima. 

Office of Telecommunications, Boulder, Colo. Inst, for 
Telecommunication Sciences. Dec 76, 60p 
PB-263 888/OPSK PC$5.25/MF$3,00 


Computers 

Windows, Harmonic Analysis and the Discrete Fourier 

Transform 

Fredric J. Harris. 

Naval Undersea Center San Diego Calif Sep 76, 70p 
AD-A034 956/3PSK PCS5.25/MF$3.00 

A Semi-Fast Fourier Transform Algorithm Over GF (2 to the 
mth power) 

Dilip V. Sarwate. 

Illinois Univ At Urbana-Champaign Coordinated Science Lab 
Sep 76, 26p 

AD-A034 982/9PSK PC$4.50/MF$3.00 


Test Procedures for Multics Security Enhancements 

M. Gasser, S. R. Ames, Jr., and L. J. Chmura. 

Mitre Corp Bedford Mass Dec 76, 331p 
AD-A034 986/OPSK PCS12,00/MF$3.00 

Software Visibility for the Program Manager 
Alan J. Driscoll. 

Defense Systems Management School Fort Belvoir Va Nov 
76, 57p 

AD-A03S 164/3PSK PCS.5.2.5/MJ-$3.00 

The Architecture of a Database Computer. Part II. The Design 
of Structure Memory and its Related Processors 
David K. Hsiao, and Krishnamurthi Kannan. 

Ohio State Univ Columbus Computer and Information 
Science Research Center Oct 76, II 3p 
AD-A035 J78/3PSK PCS6.50/MF$3.()0 

Transmission Control Protocol Specification 
Jonathan B. Postcl, Larry L. Garlick, and Raphael Rom. 
Stanford Research Inst Menlo Park Calif Augmentation 
Research Center Jul 76, 183p 
AD-A035 337/5PSK PC$9.00/MFJ3.00 

Terminal-to-Host Protocol Specification 

Jonathan B. Postel, Larry L. Garlick, and Raphel Rom. 
Stanford Research Inst Menlo Park Calif Augmentation 
Research Center Jul 76, 178p 
AD-A035 338/3PSK PC$9.00/MF$3.00 

Commercial Data Management Processor Study 
J. Cullinane, R. Goldman, T. Meurer, and R. Nawara. 
Cullinane Corp Wellesley Mass Dec 75, 84p 
AD-A035 790/5PSK PC$6.00/MF$3.U0 

Some Considerations and Models for the Distribution of a Data 
Base 

Joyce Elam, and Joel Stutz, 

Texas Univ At Austin Center for Cybernetic Studies May 76, 
32p 

AD-A035 923/2PSK PCS4.50/MF$3.00 

Software Acquisition Management Guidebook: Statement of 

Work Preparation 

J. B. Glore, and W. R. Bjerstedt. 

Mitre Corp Bedford Mass Jan 77, 104p 

AD-A035 924/OPSK PC$6.50/MF$3.00 

Data Base Management Study 

Teledyne Brown Engineering, Huntsville, Ala. Advanced 
Projects Div. May 76, 62p 
N77-15671/9PSK PCS5.25/MFS3.00 

Computer Science and Technology: Foreign and Dome.stic 
Accomplishments in Magnetic Bubble Device Technology 
Robert B. J. Warner, and Peter J. Calomeris. 

National Bureau of Standards, Washington, D.C. Inst, for 
Computer Sciences and Technology. Jan 77, 54p 
PB-263 123/2PSK PCS5.25/MFS3.00 

Accessing Individual Records from Personal Data Files Using 
Non-Unique Identifiers 

Gwendolyn B. Moore, John L. Kuhns, Jeffrey L. Trefftzs, 
and Christine A. Montgomery. 

Operating Systems, Inc,, Woodland, Calif. Feb 77, 206p 
PB-263 176/OPSK PC$9.25/MF$3.00 


Technical Profile of Seven Data Element Dictionary/Directory 
Systems 

Belkis Leong-Hong, and Beatrice Marron. 

National Bureau of Standards, Washington, D.C. Feb 77, 

49p 

PB-263 177/8PSK PCS4.50/MF$3.00 

A Contribution to Computer Typesetting Techniques, 9 track 
ASCII 

Robert C. Thompson. 

National Bureau of Standards, Washington, D.C. Office of 
Standard Reference Data. Apr 76, mag tape Source tape is 
in ASCII character set. Character set restricts preparation to 
9 track one-half inch tape only. Identify recording mode by 
specifying density only. Call NTIS Computer Products, if you 
have questions. Price includes documentation, PB-251 845. 
PB-263 925/OPSK Mag Tape $300.00 ;Foreign $480.00 


Electrotechnology 

Lightning, Surge, and Transient Protection (A Bibliography 
with Abstracts) 

William E. Reed. 

National Technical Information Service, Springfield, Va. Mar 
77, 103p 

NTIS/PS-77/0099/0PSK PC$28.00/MF$25.00 


Energy 

Site Energy Handbook. Volume I. Methodology for Energy 
Survey and Appraisal 

Energy Research and Development Administration, 
Washington, D.C. Div. of Solar Energy. Oct 76, 172p 
ERDA-76-131-1PSK PC$8.00/MF$3.00 

Preparation of a Coal Conversion Systems Technical Data 
Book. Final Report, October 31, 1974— April 30, 1976 
Institute of Gas Technology, Chicago, 111. 1976, 620p 
FE-I730-21/PSK PC$i6.25/MFS3.00 

Diesel Exhaust Emissions (Citations from the NTIS Data Base) 
Edward J. I,ehmann. 

National Technical Information Service, Springfield, Va. Mar 
77, 145p 

NTIS/PS-77/0104/8PSK PC$28 .00/MFS25 .00 

Natural Gas Deregulation Analysis 

John Neri, David Nissen, Mark Rodekohr, James Sweeney, 
and William Hogan. 

Federal Energy Administration, Washington, D.C. Jan 76, 

63p 

PB-261 599/5PSK PCS5,25/MF$3.00 

EPA Program Conference Report Fuel Cleaning Program: 

Coal 

Z. S. Altschuler, R. R. Hucko, G. A. Isaacs, P. S. Jacobson, 
and J. D. Kilgore. 

Environmental Protection Agency, Washington, D.C. Office 
of Energy, Minerals and Industry. Oct 76, 127p 

PB-262 876/6PSK PC$7.25/MF$3.00 

Methanol as an Automotive Fuel: A Summary of Research in 
the M.l.T. Energy Laboratory 

Richard G. Donnelly, John B. Heywood, Jules LoRusso, 
Frank 07Brien, and Thomas B. Reed. 



OF THE MONTH 


COBOL COMPILER VALIDATION 
SYSTEM 

Validates COBOL compilers and ensures their 
conformance to the Federal Standard prescribed in 
FIPS PUB 21-1. Consists of: 

Audit Routines — Includes tests and procedures to 
indicate test results. 

Executive Routine — Provides a method for resolv- 
ing implementor names, updating source pro- 
grams, generating operating system control cards, 
and compiling audit routines in a particular en- 
vironment. 

Order AD/A-036 173/CPP. $550 North American 
continent only. For more information on this and 
other computer products, call or write the NTIS 
Computer Products office, (703) 557-4763, 5285 Port 
Royal Road, Springfield, VA 22161. 


Massachusetts Inst, of Tech., Cambridge, Energy Lab. Apr 
76, 57p 

PB-262 980/6PSK PC$5.25/MF$3.00 

Energy Management Guide for Light Industry and Commerce. 
EPIC Energy Management Series 
W. J. Kelnhofer, and L. A. Wood. 

National Bureau of Standards, Washington, D.C. Office of 
Energy Conservation, Dec 76, 3 Ip 
PB-263 121/6PSK PC$4.50/MF$3.00 

Chemical Kinetics of the Gas Phase Combustion of Fuels (A 
Bibliography on the Rates and Mechanisms of Oxidation of 
Aliphatic Cl to CIO Hydrocarbons and of Their Oxygenated 
Derivatives 

Francis Westley. 

National Bureau of Standards, Washington, D.C. Chemical 
Kinetics Information Center. Oct 76, 142p 

PB-263 122/4PSK PC$7.25/MF.S3.0Q 

Assessment of a Single-Family Residence Solar Heating System 
in a Suburban Development Setting. Project Phoenix 
J. D. Phillips. 

Colorado Springs Dept, of Public Utilities, Colo. Jul 76, 

180p 

PB-263 192/7PSK PC$9.00/MF$3.00 

Barriers to the Use of Wind Energy Machines: The Present 
Legal/Regulatory Regime and a Preliminary Assessment of 
Some Legal/Political/Societal Problems 
Rita Falk Taubenfeld, and Howard J. Taubenfeld. 

Societal Analytics Inst., Inc., Dallas, Tex. Jul 76, 159p 

PB-263 576/lPSK PC$8.00/MF$3,00 



The Correlation of Peak Ground Acceleration Amplitude with 
Seismic Intensity and Other Physical Parameters 
L. J. O’Brien, J. R. Murphy, and J. A. Cahoud. 

Computer Sciences Corp., Falls Church, Va. Mar 77, 96p 
PB-263 972/2PSK PC $6. 00/M F$3. 00 

Cooling Tower Energy Studies - Conservation Techniques 
Applicable to Existing Installations Plus Comparative 
Economics and Energy Requirements of Mechanical and 
Natural Draft Towers 

Gordian Associates, Inc., New York. Feb 77, 70p 
PB-264 028/2PSK PC$5.25/MFS3.00 


Environment 


Synthetic Fuels from Municipal, Industrial, and Agricultural 
Wastes (Citations from the NTIS Data Base) 

Audrey S. Hundemann, and Edward J. Lehmann. 

National Technical Information Service, Springfield, Va. Mar 
77, 88p 

NTIS/PS-77/0112/1PSK PC$28.00/MF$25.00 

Dispersion Equation Solutions hy Calculator. A Guide for Air 
Pollution Engineers and Scientists 

Richard A. Porter. 

Texas Air Control Board, Austin. Nov 75, 69p 
PB-262 135/7PSK PCS5.25/MF$3.00 

The Effect of Gypsum on Lime-Fly Ash Compositions: 
Literature Review and Annotated Bibliography 

L. M. Smith, and A. Kawam. 

Gillette Research Inst., Inc., Rockville, Md. May 76, 37p 
PB-262 531/7PSK PC$4.50/MF$3.00 

Effect of a Flyash Conditioning Agent on Power Plant 

Emissions 

L. E. Sparks. 

Industrial Environmental Research Lab., Research Triangle 
Park, N.C. Particulate Technology Branch. Oct 76, 136p 
PB-262 602/6PSK PC$7.25/MF$3.00 

Proceedings: Symposium on Flue Gas Desulfurization-New 

Orleans, March 1976. Volume II 

R. D. Stern, W. II. Ponder, and R. C. Christian. 

Industrial Environmental Research Lab., Research Triangle 

Park, N.C. May 76, 453p 

PB-262 722/2PSK PC$14.50/MF$3.00 

PCB Emissions from Stationary Sources: A Theoretical Study 
Herman Knieriem, Jr. 

Monsanto Research Corp., Dayton, Ohio. Dayton Lab. Oct 
76, 43p 

PB-262 850/lPSK PC$4.50/MF$3.00 

Acute Toxicity of Selected Organic Compounds to Fathead 
Minnows 

Vincent R. Mattson, John W. Arthur, and Charles T. 
Walbridge. 

Environmental Research Lab., Duluth, Minn, Oct 76, 20p 
PB-262 897/2PSK PC$4.00/MF$3.00 

Hazardous Waste Management Facilities in the United States - 
1977 

Matthew A. Straus. 

Environmental Protection Agency, Washington, D.C. Office 


of Solid Waste Management Programs. Jan 77, 66p 
PB-262 917/8PSK PC$5.25/MF$3.00 

Evaluation of the General Motors’ Double Alkali S02 Control 
System 

Edward Interess. 

Little (Arthur D.), Inc., Cambridge, Mass. Jan 77, 99p 
PB-263 469/9PSK PCS6.00/MF$3,00 

Tall Stacks and the Atmospheric Environment 

Jeremy M. Hales. 

Battelle Pacific Northwest Labs., Richland, Wash. Mar 76, 
178p 

PB-263 561/3PSK PC$9.00/MF$3.00 

Procedures for Cascade Impactor Calibration and Operation in 
Process Streams 
D. Bruce Harris. 

Industrial Environmental Research Lab., Research Triangle 

Park, N.C. Jan 77, 123p 

PB-263 623/lPSK PC$6.50/MF$3.00 

Public Comments and Task Force Responses Regarding the 
Environmental Survey of the Reprocessing and Waste 
Management Portions of the LWR Fuel Cycle. (NUREG-0116) 
Nuclear Regulatory Commission, Washington, D.C. Office of 
Nuclear Material Safety and .Safeguards. Mar 77, 427p 
PB-264 035/7PSK PCS 14.00/MF$3.00 


Medical Screening Programs 
for Infants 

The proceedings of a 1972 symposium on laboratory 
screening techniques for newborn and selected 
high-risk infants are now available from NTIS. The 
symposium was jointly sponsored by the U.S. and 
Poland and was attended by representatives of 23 
countries — scientists who developed original screen- 
ing techniques, people from diagnostic centers 
working on automation of screening methods, and 
experienced practitioners. Participants shared the 
findings of new or modified techniques, hoping to 
improve the effectiveness and stimulate the growth 
of laboratory and mass screening programs for new- 
born infants. Subject matter addressed covered 
laboratory screening techniques applied to dried 
blood spots, capillary blood, urine samples, and 
other biological specimens. Participants also dis- 
cussed the organization and administrative aspects 
of laboratory and mass screening programs. 
Proceedings of the International Symposium on 
Laboratory Screening Techniques for Inborn Errors 
of Metabolism in Newborn and Selected High-Risk 
Infants lists the names and addresses of all sym- 
posium participants and offers a convenient subject 
index. Order the Proceedings, PB-263 027-T/PSK, at 
$11.00 per paper copy and $3.00 per microfiche copy 
from the National Technical Information Service, 
5285 Port Royal Road, Springfield, Va. 22161. 
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Order Here 
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National Technical Information Service 
U.S. DEPARTMENT OF COMMERCE 
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NTIS products and services, please fill out the In- 
formation Requested section in full because this 
section is sent to a separate office to fill your in- 
formation requests. 

To complete the Items Requested section first 
enter your name and address on the space 
provided. Fill out the Ship to section only if you 
wish your order shipped to another address than 
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When you indicate the method of payment, 
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panied by payment, you will be billed an ad- 
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continent; $3,00 outside North American charge 
per item. Just check the Priority mail box, If 
you’re really pre,ssed for time, call the NTIS 
Rush Flandling Service (703) 557-4700. For a 
$10,00 charge per item, your order will be sent 
priority mail within 24 hours. Or, you can pick 
up your order in the Washington Information 
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Health Planning 


Reducing Excess Hospital Capacity 
Walter McClure. 

InterStudy, Excelsior, Minn. Oct 76, 236p 
HRP-0015199/3PSK PCS9.50/MF$3.00 

Library and 
Information Sciences 

Development of a Model Library Microform Center 
Philip J. Gary. 

Gary (Philip J.) and Associates, Marina del Rcy, Calif. 
Micrographic Systems Consulting Services. Aug 76, 56p 
PB-263 586/OPSK PCS5.25/MF$3.00 


Materials Sciences 

Detection of Tightly Closed Flaws by Nondestructive Testing 
(NDT) Methods in Steel and Titanium 

W. D. Rummel, R. A. Rathke, P. H. Todd, Jr., T. L. Tedrow, 
and S. J. Mullen. 

Martin Marietta Aerospace, Denver, Colo. Sep 76, 109p 
N77-1S133/0PSK PC$6..^0/MF$3.00 



OF THE MONTH 


New Improved Method for Synthesis 
of Methotrexate 

Consists of three process improvements in the 
production of an anticancer agent, antifolate metho- 
trexate — 

• Selection of a PH value that minimizes the for- 
mulation of pteridines 

• Elimination of the sodium acetate or acetate 
buffer previously used 

• Conversion of the sulfite acid addition salt to 
hydrochloride salt 

Results ... an easier process at about half the 
current cost . . . with better yields and purer prod- 
ucts. 

For more information on this and other Govern- 
ment-owned technology, contact George Kudravetz 
at (703) 557-4733 or write him at NTIS, 5285 Port 
Royal Road, Springfield, VA 22161. 


Encapsulation by Entrapment 

Baruch S. Shasha, William M. Doane, and Charles R. Russell. 
Department of Agriculture, Washington, D.C. Filed 19 Oct 
76, 42p This Government-owned invention available for 
U.S. licensing and, possibly, for foreign licensing. Copy of 
application available NTIS. 

PAT-APPL-733 968/PSK PC$4.50/MF$3.00 

The Electroslag Melting Process 
R. H. Nafziger. 

Bureau of Mines, Albany, Oreg. Albany Metallurgy Research 

Center. 1976, 234p 

PB-263 219/8PSK PC$9.50/MF$3.00 

Medicine and Biology 

Technical Notes: MEDLARS Indexing Instructions: 
Supplements V. Enzyme Key and VI. Fungu-S Key 
Elizabeth J. VanLcntcn. 

National Library of Medicine, Bethesda, Md. Feb 77, 34p 
PB-263 335/2PSK PC$4.50/MF$3.00 

Ocean Technology 
and Engineering 

standardized Stern Project. Volume I 

Charles H. Gross, Jr. 

Henry (J. J.) Co., Inc., New York. Sep 76, 59p 
PB-262 926/9PSK PC$5.25/MF$3.00 

Standardized Stern Project. Volume II. Executive Summary 
Charles H. Gross. Jr. 

Henry (J. J.) Co., Inc., New York. Sep 76, lOp 
PB-262 927/7PSK PC.S4.00/MF$3.00 


T ransportation 

Transportation System Safety Methodology 

M. Horodniceanu, E. J. Cantilii, M. Shooman, and E. J. 

Pignataro. 

Polytechnic Inst, of New York, Brooklyn, Nov 76, 147p 
PB-262 793/3PSK PC$7.2.VMF$3,00 

Impact of a Suburban Rapid Transit Line of Fuel Consumption 
and Cost for the Journey-To-Work. Analysis of the 
Philadciphia-Lindenwold High-Speed Line 

David E. Boyce, Khanh Nguyen, Thierry Noyelle, and Kent 
Webb. 

Wharton School of Finance and Commerce, Philadelphia, Pa. 
Dept, of Regional Science. Dec 75, 74p 
PB-263 048/lPSK PC$5.25/MFS3.00 

Carpool Incentives: Evaluation of Operational Experience 

Cambridge Systematics, Inc., Mass. Mar 76, 176p 
PB-263 050/7PSK PC$9.00/MF$3.00 

Aesthetic and Environmental Considerations in the Design of 
Elevated Transportation Structures 

M. L. Silver, and L. S. Bell. 

Illinois Univ. at Chicago Circle. Dept, of Materials 

Engineering. Jun 75, 134p 

PB-263 158/8PSK PC$7.25/MF$3.00 

Federal Support for the Development of Alternative 
Automotive Power Systems: The General Issue and the 
Stirling, Diesel, and Electric Cases 

l.awrence H. I.inden, John B. Heywood, H. D. Jacoby, and 
Howard Margolis. 

Massachusetts Inst, of Tech., Cambridge. Energy Lab. Mar 
76, 406p 

PB-263 523/3PSK PC$13,25/MF$3.00 
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Dear Friend, 

Are you looking for a great gift idea or a special 
addition to your own library? Send now for the 
History of the Star Spangled Banner from 1814 to 
the Present. Young and old will enjoy more than 
500 pages of fascinating details on the anthem’s 
publication as a poem, adaptation to music, the 
development of its popularity, inspirational meaning, 
and more. You can now purchase this attractive, 
hard cover report written in the spirit of our nation’s 
heritage for just $6.75. 

NTIS has slashed the price from $13,50 to $6.75 for 
a limited time only. Don’t let this opportunity pass 
you by. Act now. Order PB-196 294/PSK for $6.75. 

Ride the Subway to NTIS 

Now that the Washington Metrorail System has 
been expanded, NTIS customers will benefit from 
more convenient access to our Downtown Book- 
store. Located just two short blocks from Metro 
Center — the heart of the subway system— the Book- 
store is just minutes away from many points in and 
around the city. 

Whether you’re an out-of-towner traveling the 
Blue Line from National Airport or Union Station 
... or a local customer riding the Red Line across 


town, be sure to visit the NTIS Bookstore for those 
last-minute orders or rush pick-ups. 

Use Code-A-Phone After 5 

Another new customer convenience recently 
added by NTIS is code-a-phone . . . enabling you to 
place orders, leave messages, or ask questions even 
after our regular Bookstore hours are over. The 
very next day of business, an NTIS representative 
will return your call with news about your order or 
with an answer to your question. Keep the code-a- 
phone number — (202) 724-3383 — in a handy place 
at your home and office, and the next time you need 
something, give us a call. 

Sincerely, 


Rise L. Burggraff 
Editor 

P.S. The Premium Service toll-free number was 
incorrectly announced in an earlier catalog. Please 
note the following numbers for Premium Service 
orders. (800) 325-6000; in Missouri, (800) 342-6700. 


Parkinson^s Disease and Related 
Disorders — 

Worldwide References Available 


The National Institute of Neurological and Com- 
municative Disorders and Strokes publication, 
Parkinson's Disease and Related Disorders: Citations 
from the Literature, is now available by annual sub- 
scription from NTIS. Issued monthly, the 
bibliography provides current, selected, 
systematically categorized citations from this 
sizeable body of worldwide biomedical literature. 


The National Library of Medicine's MEDLARS 
program is used for selecting appropriate items. 
Citations to "core" documents — those primarily 
clinical in approach and those connecting research 
directly to a specific disorder are included. These 
detai I Parkinson's d isease, Huntington's and 
Sydenham's choreas, dystonia musculorum 
deformans, hepatolenticular degeneration, multiple 
tics, and torticollis. 

An annual subscription to Parkinson's Disease and 
Related Disorders, costs $30. To help keep your 
specialized medical knowledge current, send your 
order to National Technical Information Service, 
5285 Port Royal Rd., Springfield, Va. 22161. 
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